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Abstract 


Nanoscale  intergranular  “glassy”  films  often  control  the  fabrication  and  mechanical 
properties  of  high-temperature  structural  ceramics.  This  AFOSR  Young  Investigator  program 
unequivocally  demonstrated  the  high-temperature  formation  of  analogous  liquid-like  grain 
boundary  (GB)  films  in  metallic  refractory  alloys.  The  bulk  computational  thennodynamic 
(CalPhaD)  methods  were  extended  to  GBs,  predicting  the  onset  of  GB  disordering  at  as  low  as 
60-85%  of  the  bulk  solidus  line.  Combined  experimental  and  modeling  studies  of  both  W  and 
Mo  based  systems  unambiguously  demonstrated  that  the  mysterious  solid-state  activated 
sintering  is  due  to  the  increased  mass  transport  in  impurity-based  liquid-like  GB  films  that  are 
thennodynamic  stabilized  below  the  bulk  solidus  line.  Determining  this  solid-state  activated 
sintering  mechanism  solved  an  outstanding  scientific  problem  that  has  puzzled  the  materials 
community  for  >50  years.  Subsequently,  Ni-doped  Mo  was  selected  for  systematical  evaluation 
via  characterizing  well-quenched  specimens  and  thermodynamic  modeling.  Finally,  “GB 
‘phase’  diagrams”  were  developed  as  a  new  materials  science  tool  to  control  microstructural 
evolution  and  forecast  high-temperature  properties;  their  correctness  has  been  quantitatively 
validated  by  measuring  GB  diffusivities,  direct  HRTEM  and  Auger  characterizations,  and 
atomistic  simulations.  The  developed  high-temperature  interfacial  thermodynamic  models  can 
be  further  extended  to  more  complex  multicomponent  metals  and  ceramics  in  future  studies. 
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1.  Overview 


Nanoscale,  impurity-based,  intergranular  films  (IGFs)  are  commonly  observed  in  S13N4,  SiC, 
SiAlON,  AI2O3,  ZrC>2,  BN,  B4C  and  other  high-temperature  structural  ceramics.  They  often 
control  sintering,  grain  growth,  and  mechanical  properties  [1].  With  the  support  of  this  AFOSR 
Young  Investigator  Program  (AFOSR- YIP),  we  discovered  and  studied  an  analogous  interfacial 
phenomenon  in  binary  refractory  metals: 

•  We  demonstrated  the  stabilization  of  subsolidus  quasi-liquid  IGFs  (Fig.  1)  [2]. 

>  Despite  various  indications  [1],  this  is  the  first  direct  high-resolution  transmission 
electron  microscopy  (HRTEM)  observation  of  quasi-liquid  IGFs  in  metallic  alloys. 

>  This  observation  deepens  the  understandings  of  this  class  of  interfacial  phenomena. 

•  We  demonstrated  that  short-circuit  diffusion  in  subsolidus  quasi-liquid  IGFs  result  in  “solid- 
state  activated  sintering”  in  refractory  metals.  This  work  solved  an  outstanding  scientific 
problem  that  has  puzzled  the  materials  community  for  over  50  years  [2,  3]. 

•  We  extended  bulk  computational  thermodynamic  (CalPhaD)  methods  to  model  grain 
boundaries  (GBs)  and  predicted  that  the  nanoscale  quasi-liquid  IGFs  can  be  stabilized  at 
temperatures  as  low  as  60-85%  of  the  bulk  solidus  lines  [3,  4]. 

>  Such  quasi-liquid  IGFs  can  critically  affect  the  microstructural  evolution  and  high- 
temperature  properties  [3]. 


Well-Known 


Our  Discovery 


Scientific  Significance:  Metallic  counterparts 

•  Relative  simple  systems  high-T interfacial  thermodynamics 


Help  to  understand  more  complex  ceramics,  e.g.,  Si3N4  &  UHTM  j\ 


Technological  Importance 

•  GB  diffusivity  t  (sintering?  Coble  creep?) 

•  GB  mobility?  (grain  growth?  microstructural  evolution?) 

•  High-T  mechanical  properties,  e.g.  GB  sliding?  creep? 

•  Oxidation  or  corrosion  resistance? 

•  GB  embrittlement?  Liquid-metal  embrittlement? 


Relevant  to  high- 
T  materials 
«-  (because  its 
stability  scales 
with  TJ 


A  Long-Range  Scientific  Goal:  GB  “phase”  diagrams 

♦  A  new  tool  for  mechanism-informed  design  of  high-T  materials. 


Fig- 1  :  The  well-known  equilibrium-thickness  IGFs  in  ceramic  materials  and  analogous  quasi-liquid  IGFs 

in  refractory  metals. 
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We  further  proposed  a  long-range  scientific  goal  of  developing  “GB  diagrams”  as  a  new 
tool  for  mechanism-informed  materials  design.  Fig.  2  represents  two  of  such  GB  diagrams 
developed  and  verified  in  this  AFOSR-YIP.  These  GB  diagrams  can  be  used  to  a)  provide 
useful  information  for  designing  optimal  fabrication  protocols  and  b)  forecast  certain  high- 
temperature  material  properties. 

This  interfacial  phenomenon  is  more  important  for  high-temperature  materials,  such  as 
refractory  metals  and  ceramics.  This  is  because  GB  energy  increases  with  increasing  melting 
temperature,  and  greater  GB  energies  can  stabilize  liquid-like  GB  structures  over  larger 
undercooling  ranges  and  critically  affect  microstructural  evolution  (e.g.,  sintering  and  grain 
growth)  and  high-temperature  properties  (e.g.,  creep).  Furthermore,  impurity-based  IGFs  will 
often  be  retained  upon  cooling  and  affect  low-temperature  properties  (e.g.,  embrittlement)  [1]. 


Verified  by  HRTEM: 
hE Q  =  £-ln (A/£) »  0.6  nm 
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Fig.  2  :  “GB  diagrams”  for  (a)  Ni-doped  W  and  (b)  Ni-doped  Mo,  in  which  computed  thicknesses  of  stable  quasi¬ 
liquid  IGF  (2,’s;  red  dashed  lines)  are  plotted  in  the  binary  bulk  phase  diagrams.  The  colors  represent  the 
thermodynamic  tendency  for  GBs  to  disorder.  These  computed  GB  diagrams  were  validated  via  experiments  (by 
FIRTEM,  Auger  and  GB  diffusivity  measurements)  and  atomistic  simulations  (by  U.  Alberta  collaborators). 
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2.  Background 


2.1.  Premelting 

Premelting,  first  postulated  by  Faraday  in  1842,  refers  to  the  stabilization  of  nanoscale  liquid¬ 
like  interfacial  layers  below  the  bulk  melting  temperature.  Surface  premelting  has  been  observed 

in  many  unary  solids  [18-20].  For  example,  ice  premelting  (Fig.  3)  plays  important  roles  in  snow 

* 

“sintering”,  glacier  motion  (creep)  [20],  and  the  ice-skating  process  [21]  . 

Materials  researchers  have  also  sought  to  confirm  the  existence  of  GB  premelting.  In  1989, 
Balluffi  et  al.  reported  [22]  that  GB  premelting  did  not  occur  until  0.999 rnie|ling  for  pure  Al.  The 
occurrence  of  GB  premelting  in  colloidal  crystals  was  reported  in  2005  (Fig.  4)  [23]. 

Nonetheless,  the  importance  of  GB  premelting  in  unary  materials  remains  controversial. 


In  multicomponent  materials,  however,  a  class  of  "premelting-like"  films  can  be  stabilized  at 
GBs  over  wider  ranges  of  undercooling,  wherein  GB  disordering  is  enhanced  by  a  concurrent 
adsorption  in  a  "prewetting"  regime.  Here,  “ prewetting ”  refers  to  an  adsorption  transition 
initially  proposed  in  the  Cahn  critical-point-wetting  model  [24]  for  binary  de-mixed  liquids. 


Fig.  3:  Surface  premelting  in  ice  (atomistic  simulation  by  Furukawa  and 

Nada,  1997). 


Fig.  4:  GB  premelting  in  a 
colloidal  crystal,  an  image  featured 
on  the  cover  of  Science  [23]. 


2.2.  Equilibrium-Thickness  Intergranular  Films  (IGFs) 

Researchers  studying  ceramics  have  discovered  the  widespread  existence  of  a  unique  class  of 
intergranular  films  (IGFs)  that  exhibit  the  following  distinct  characteristics  [1]: 

•  self-selecting  or  “equilibrium”  thickness  on  the  order  of  1  nm; 

•  composition  that  is  neither  observed,  nor  stable,  as  a  bulk  liquid/glass  phase  (e.g.,  the 
average  film  composition  can  lie  within  a  bulk  immiscible  gap!);  and 

•  quasi-liquid  structure  that  is  neither  completely  crystalline,  nor  fully  amorphous. 

These  nanoscale  IGFs  can  be  alternatively  understood  to  be: 

•  liquid-like  interfacial  films  that  adopt  an  equilibrium  thickness  in  response  to  several 
attractive  and  repulsive  interfacial  forces  (the  Clarke  model)  [25-27]  or 

•  a  unique  class  of  disordered  multilayer  adsorbates  (the  Cannon  model)  [28]. 

The  stabilization  of  impurity-based,  quasi-liquid  IGFs  at  subsolidus  temperatures  have  been 
observed  [2,  29-31],  provoking  an  analogy  to  “premelting”  in  unary  systems. 


The  popular  "pressure  melting"  theory  for  ice-skating  is  flawed  [21].  An  ice  skater  can  exert  a  pressure  to  cause 
local  melting  at  only  a  few  degrees  below  zero.  Premelting  and  “frictional  heating”  enable  skating  at  -35  °C. 
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In  2006,  Tang,  Carter,  and  Cannon  (MIT  and  Berkeley)  [1]  proposed  to  explain  subsolidus 
IGFs  in  binary  systems  from  coupled  GB  premelting  (structural  disordering)  and  prewetting 
(adsorption)  transitions  [32].  The  Tang-Carter-Cannon  model  [32]  serves  as  a  qualitative  basis 
for  understanding  IGFs  in  ceramics.  Electrostatic  [26]  and  van  der  Waals  (vdW)  London 
dispersion  forces  [33]  must  be  separately  included,  resulting  in  complex  adsorption  and  wetting 
behaviors.  For  example,  nanoscale  equilibrium-thickness  IGFs  can  often  persist  above  the  bulk 
solidus  line,  due  to  the  presence  of  the  attractive  van  der  Waals  (vdW)  London  dispersion  forces 
in  ceramics  [1,  34].  IGFs  in  ceramics  can  also  represent  a  metastable  equilibration. 

In  addition  to  the  well-known  silicate-based  IGFs  in  Si3N4,  SiC  and  other  structural  ceramics, 
thin  interfacial  films  of  similar  characteristics  have  been  found  in  an  increasing  number  of  other 
material  systems  (Fig.  5)  [1],  including: 

•  non-silicate-based  IGFs  in  various  oxide  ceramics; 

•  IGFs  at  ceramic  hetero-interfaces; 

•  silicate-based  IGFs  at  metal-oxide  interfaces; 

•  free-surface  counterparts  to  IGFs;  and 

•  metallic  IGFs. 

The  PI  has  recently  published  a  43-page  Critical  Review  article  about  the  observations, 
theories,  and  technological  importance  of  nanoscale  quasi-liquid  interfacial  films  [1].  More 
specifically,  the  existence  of  metallic  IGFs  was  unequivocally  and  directly  confirmed  in  our 
AFOSR-YIP  [2-4,  9].  Moreover,  understanding  “surface  amorphous  films”  (SAFs,  i.e.,  the  free- 
surface  counterparts  to  IGFs)  [13]  can  help  to  control  nanoscale  native  oxide  films  and  their 
detrimental  effects  in  the  densification  of  non-oxide  high-temperature  materials. 


Hetero-Interfaces 
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Free  Surfaces 


Oxide -Oxide 
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Ackler  &  Chiang, 


film 
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Equilibrium-Thickness:  Clarke,  1987 


Metallic  IGFs 


Luo  et  al.,  2005 


Pure  Ice  Surface: 

Faraday,  1860's 

Enhanced  by  Minor 
Impurities,  e.g.  NaCI 

(Wettlaufer,  1999) 

Surface  ->  GB  & 
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Widespread  Existence 
of  Similar  IGFs: 


Interface 

(Avishai  et  al,  2005) 


Kwonefa/.,  1999 
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SiC,  AIN,  SiAlON, 
AljOj,  BaTi03, 
SrTi03,  PbjRujO,, 
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or  impurities) 


No  Si02,  T<  Teuteotic 
Wang  &  Chiang,  1998 
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Fig.  5:  The  widespread  existence  of  various  impurity-based  quasi-liquid  interfacial  fdms.  Reprinted  from  the  Pi's 

critical  review  (an  AFOSR-supported  publication  [1]). 
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3.  Summary  of  Key  Results 

The  experimental  and  computational  methods  used  in  this  project  are  schematically 
illustrated  in  Fig.  6.  The  detailed  experimental  procedures  are  documented  in  Appendix  A.  Key 
results  are  discussed  in  the  following  sections. 


Methods 


Sample 

bolder 


Experiments 

•  Carefully-prepared  specimens 

•  High-purity 

♦  Isothermally  annealed 


S 


A  barrrl  of 
eolil'Mutcr 


..Sample 
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body 


„  Quench 
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TT 

I 


Quench  oil 


•  Well  quenched  (using  custom-built  furnace) 
Characterized  by  HRTEM,  Auger  etc. 

Controlled  sintering  experiments  A  GB  diffusivity 


Quench  time  <  Is 


Theories  &  Modeling 


•  Extend  bulk  CalPhaD  (Calculation  of  Phase  Diagram)  to  GBs 

•  Statistical  thermodynamic  models  A  interfacial  energies 

•  Atomistic  modeling  (via  collaboration  w /  U.  Alberta) 


Fig.  6:  Summary  of  experimental  and  theoretical  methods  used  in  this  study. 


3.1.  Nanoscale  Quasi-Liquid  IGFs  in  Doped  W 

We  unequivocally  demonstrated  the  stabilization  of  nanoscale  quasi-liquid  IGFs  well  below 
the  bulk  eutectic  temperature  in  Ni-doped  W  (Fig.  7)  [2],  While  indirect  evidence  has  been 
reported  previously  [42],  this  is  the  first  direct  HRTEM  evidence  supporting  the  existence  of 
such  nanoscale  quasi-liquid  IGFs  in  metallic  alloys. 


Fig.  7:  Representative  HRTEM  image  of  nanoscale  quasi-liquid  IGFs  that  are  thermodynamically  stable  at  95  K 

below  the  bulk  eutectic  temperature  in  Ni-doped  W. 

Reprinted  from  an  AFOSR  supported  publication  in  Acta  Materialia  [2], 
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By  extending  bulk  CalPhaD  (Calculation  of  Phase  Diagram)  methods  to  GBs,  we  predicted 
that  quasi-liquid  IGFs  can  be  stabilized  at  W  GBs  as  low  as  60-85%  of  the  bulk  solidus 
temperatures  [4],  The  basic  concepts  of  this  model  are  schematically  illustrated  in  Fig.  8.  A 
refined  thennodynamic  model  is  elaborated  in  Appendix  B.  The  key  results  of  this  model  [3,  4] 
are  briefly  discussed  below. 


Thermodynamic  Principle 

Crystal 


I  subsolidus  q 
thermodynamically  stable  if: 


forming 

undercooled  liquid  ^ 
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Combining 
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Interfacial  Interfacial 
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/(«,)  =  ! 


Equilibrium  thickness  (hEQ): 
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dh 


A  Thermodynamic  Variable:  X 


(a)  (b) 

Fig.  8  :  Schematic  illustration  of  the  proposed  thermodynamic  model. 


As  shown  in  Fig.  8(a),  a  subsolidus  quasi-liquid  IGF  can  be  thermodynamically  stable  if: 

Mmwi/h<(-Ar)=/*-2rd  (i) 

Where:  h  is  the  film  thickness,  ygb(0)  is  the  excess  energy  of  a  random  GB  without  adsorption,  and 
yc i  is  the  crystal-liquid  interfacial  energy.  The  volumetric  free  energy  to  form  an  undercooled 
liquid  (AG’amorph.)  of  a  binary  alloy  can  be  quantified  using  the  bulk  CalPhaD  methods.  A 
thermodynamic  variable  A  is  defined  to  represent  the  maximum  thickness  of  a  quasi-liquid  IGF 
that  can  be  thermodynamically  stable  (without  consideration  of  interfacial  forces): 

/l  =  Max{-Ay /  AGamorph}  ,  (2) 

The  definition  and  quantification  of  the  thermodynamic  variable  X  are  schematically  illustrated 
in  Fig.  8(b).  Further  details  of  this  and  refined  models  are  discussed  in  Appendix  B. 

Five  computed  “GB  diagrams”  are  shown  in  Fig.  2  (Ni-doped  W;  Ni-doped  Mo)  and  Fig.  9 
(Pd-doped  W;  Co-doped  W;  and  Fe-doped  W),  where  the  lines  of  constant  X  are  plotted  in  the 
bulk  binary  phase  diagrams. 

For  metallic  alloys,  the  dominance  of  an  exponentially-decaying  interfacial  force  allows  a 
simple  estimate  of  the  actual  “equilibrium”  IGF  thickness  as:  h\_Q  ~  c-ln(X/c),  where  is  the 
coherent  length.  The  “GB  solidus  temperature”  (7gbs)  can  then  be  estimated  by  /1(7gbs)  = 

0.5-1  nm.  Interfacial  forces  are  more  complex  in  ceramic  materials  [1]. 

In  metallic  alloys,  the  film  thickness  is  divergent  as  the  bulk  solidus  line  is  approached  from 
below,  above  which  a  bulk  liquid  phase  appears  and  completely  wets  the  GBs.  In  ceramics, 
however,  attractive  vdW  forces  can  stabilize  nanoscale  IGFs  above  the  bulk  solidus  lines. 
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Fig.  9:  Computed  GB  diagrams  for  (a)  Pd-doped  W,  (b)  Co-doped  W,  and  (c)  Fe-doped  W.  Two  additional 
computed  GB  diagrams  for  Ni-doped  W  and  Ni-doped  Mo  are  shown  in  Fig.  2.  In  each  bulk  phase  diagram,  a  series 
of  dashed  and  dotted  lines  were  plotted  to  represent  the  conditions  that  4,  2,  1  and  0.5  nm  thick  quasi-liquid  IGFs 
can  be  thermodynamically  stable  at  GBs.  The  ranges  of  subsolidus  activated  sintering  are  also  labeled. 
Redrawn  from  an  AFOSR  supported  publication  in  Curr.  Opin.  Solid  State  Mater.  Sci.  [3]. 


3.2.  Solid-State  Activated  Sintering:  Resolving  a  50-Year-Old  Mystery 

Solid-state  (subsolidus)  activated  sintering  refers  to  the  enhancement  of  densification  due  to 
solid-state  sintering  additives.  In  their  classic  1978  review  [43],  Coble  and  Cannon  stated  that: 

" The  most  significant  changes  which  have  been  taken  place  in  recent  years  regard 
the  finding  with  respect  to  densification  below  the  eutectics  in  numerous  systems, 
notably  WC:Co,  W:Ni,  ..." 

The  exact  mechanism  has  puzzled  the  materials  community  for  decades. 

A  prior  study  by  the  PI  provided  the  first  insight  towards  resolving  this  mystery  [44].  In 
BFCfi-doped  ZnO,  nanoscale  “glassy”  IGFs  were  found  to  be  stabilized  at  GBs  well  below  bulk 
Teutectic ,  occurring  concurrently  with  the  onset  of  activated  sintering,  implying  that  the  enhanced 
sintering  is  likely  due  to  the  short-circuit  diffusion  in  these  subsolidus  IGFs  [44],  While  a 
ceramic  is  not  the  best  system  to  unequivocally  prove  the  hypothesis,  due  to  the  complex 
interfacial  forces,  we  proved  this  hypothesized  activated  sintering  mechanism  systematically 
using  the  binary  refractory  metals  in  this  AFOSR-YIP,  as  elaborated  subsequently. 

It  has  been  commonly  understood  since  the  1950’s  that  the  addition  of  <  1  %  of  various 
transition  metals  can  significantly  increase  the  densification  of  W  and  Mo  in  the  subsolidus 
region  [43,  45-51].  In  a  classical  model,  the  solid-state  "activator"  is  presumed  to  be  the 
secondary  bulk  crystalline  phase  that  completely  wets  the  GBs.  We  recently  unequivocally 
illustrated  that  solid-state  wetting  does  not  occur  (Fig.  10;  Fig.  13)  [2,  3,  9].  Instead,  quasi-liquid 
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IGFs  can  be  stabilized  at  GBs  well  below  the  bulk  solidus  line  and  lead  to  an  enhanced  GB 
diffusion  and  activated  sintering  [2-4,  9]. 


Activated  Sintering  in  W-Ni 

Gupta,  Yoon,  Meyer  &  Luo,  Acta  Mater.  55:3131  (2007) 


Luo,Cr/t.  Rev.  Solid  State  Mater. 
Sc/.  32:67-101(2007) 


A  classical  wisdom: 

Solid-state  wetting:  the  crystalline  secondary 
phase  completely  wets  GBs 


Fig.  10:  In  contrast  to  the  classical  wisdom,  the  Ni-rich  crystalline  phase  does  not  wet  the  W  GB.  HRTEM  and 
Auger  electron  spectroscopy  (AES)  revealed  the  stabilization  of  nanoscale  Ni-enriched,  liquid-like  IGFs  below  the 
bulk  eutectic  temperature.  Combined  experimental  and  modeling  studies  unambiguously  demonstrated  that  the 
mysterious  solid-state  activated  sintering  is  due  to  the  increased  mass  transport  in  impurity-based  liquid-like  GB 
films  that  are  thermodynamic  stabilized  below  the  bulk  solidus  line.  Determining  this  mechanism  for  solid-state 
activated  sintering  in  refractory  alloys  solved  an  outstanding  scientific  problem  that  has  puzzled  the  materials 

community  for  >50  years. 


Subsolidus  activated  sintering  experiments  have  been  conducted  for  W  using  various  dopants 
with  significantly  different  effectiveness  [48,  52],  enabling  a  critical  test  of  the  hypothesized 
sintering  mechanism.  On  the  other  hand,  by  using  onset  sintering  as  an  indicator  for  GB 
disordering,  we  can  test  whether  our  model  can  predict  GB  disordering,  which  has  broader 
applications  beyond  sintering  [4],  The  comparisons  between  the  predictions  and  experiments  for 
five  systems  (W-M,  M  =  Pd,  Ni,  Co,  Fe  and  Cu)  are  illustrated  in  Table  1  and  Fig  11  [4].  These 
results  not  only  support  the  hypothesized  subsolidus  activated  sintering  mechanism,  but  also  the 
correctness  and  usefulness  of  our  model  described  in  §3. 1 . 

Table  1 :  Measured  onset  sintering  temperatures  vs.  predicted  GB  disordering  temperatures.  Reprinted  from  an 

AFOSR  supported  publication  in  Appl.  Phys.  Lett.  [4]. 


T el  p 

(K) 

^ Ye/p 
(Mm1) 

Cinter  (K) 

7gbs  (K) 

W-Pd 

2088 

-0.60 

1090  ±23 

<  1 1 4 1 a 

W-Ni 

1768 

-0.52 

1 1 50  ±  18 

1121-1470 

W-Co 

1962 

-0.40 

1301  ±49 

1140-1644 

W-Fe 

1910 

-0.37 

1308  ±50 

1273-1664 

W-Cu 

1357 

+0.49 

No  activated  sintering 

No  IGF 
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(a) 


(b) 


Fig.  11 :  (a)  Computed  IGF  thickness  vs.  temperature  for  dopant-saturated  W  specimens,  (b)  The  experimental 
sintering  (densification)  rates.  Reprinted  from  an  AFOSR  supported  publication  in  Appl.  Phys.  Lett.  [4]. 


3.3.  “GB  Diagrams”  -  A  New  Materials  Science  Tool? 

In  this  AFOSR- YIP,  we  developed  a  long-range  scientific  goal  of  establishing  “GB 
diagrams”  as  a  new  tool  for  mechanism-informed  materials  design,  which  involves  three 
concepts: 

•  Guided  by  GB  diagrams,  fabrication  protocols  can  be  designed  to  utilize  the  most 
appropriate  GB  structures  to  achieve  optimal  microstructures,  e.g.,  utilizing  liquid-like  GBs 
for  low-temperature  sintering; 

•  Heat  treatment  protocols  can  be  devised  to  adjust  GB  structures  to  achieve  the  desired 
properties,  e.g.,  to  de-wet  impurity-based  IGFs  to  mitigate  GB  embrittlement;  and 

•  GB  diagrams  provide  information  for  predicting  high-temperature  properties,  e.g.,  creep 
resistance. 

The  necessity  of  developing  such  GB  diagrams  is  demonstrated  by  our  studies  of  the  solid-state 
activated  sintering  of  ceramics  [44]  and  metals  [2],  Since  nanoscale  quasi-liquid  IGFs  can  fonn 
well  below  the  bulk  solidus  lines  and  result  in  enhanced  sintering  behaviors  similar  to  liquid- 
phase  sintering  [4],  bulk  phase  diagrams  are  not  adequate  for  designing  activated  sintering 
protocols.  On  the  other  hand,  our  results  indeed  showed  that  the  onset  activated  sintering  can  be 
predicted  from  the  computed  GB  diagrams  (Table  1;  Fig.  9;  Fig.  1 1). 

Furthermore,  such  GB  diagrams  can  lead  to  applications  in  the  mechanism-informed  design 
of  fabrication  protocols  (beyond  sintering)  and  the  prediction  of  high-temperature  properties. 


3.4.  A  Systematical  Study  of  Ni-Doped  Mo 

We  recently  systematically  evaluated  the  Ni-doped  Mo  system  via  characterizing  well- 
quenched  specimens,  in  conjunction  with  thermodynamic  modeling.  In  contrast  to  prior  beliefs 
[50,  53],  we  found  that  the  S-NiMo  phase  does  not  wet  Mo  GBs  in  the  solid  state  (Fig.  13).  On 
the  other  hand,  the  Ni-rich  liquid  wets  Mo  GBs  completely  above  the  solidus  line  (Fig.  13). 
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Furthermore,  HRTEM  and  AES  analysis  demonstrated  that  nanometer-thick  quasi-liquid  films 
persist  at  GBs  into  the  single-phase  region  (X <  Xsoiidus),  where  the  bulk  liquid  phase  is  no  longer 
stable  (Fig.  12).  While  such  quasi-liquid  IGFs  promote  sintering  at  firing  temperatures,  they  will 
become  brittle  5-NiMo  layers  upon  cooling,  resulting  in  GB  embrittlement  [50]. 

Furthennore,  molecular  dynamic  simulation  conducted  by  our  collaborator,  Professor  Hao 
Zhang  from  the  University  of  Alberta  in  Canada,  supported  our  experiments  in  proving  that  the 
addition  of  Ni  (below  the  bulk  solid  solubility  limit)  causes  GBs  to  widen  and  become  more 
disordered  (Fig.  14).  This  also  significantly  increases  the  GB  diffusion  (Fig.  14(c)). 


Measured:  h  =  0.80  ±  0.12  nm 


Mo  +  1  at.  %  Ni.  Quenched  from  1495  °C 

(a) 


30  sec  Sputter  Time 
=  1  nm  for  a  Si02 
reference  specimen 


E 

o 

o 


(b)  (c) 

Fig.  12:  (a)  HRTEM  images  and  (b)  AES  analysis  revealed  the  thermodynamic  stabilization  of  ~  0.8  nm  thick 
quasi-liquid  IGFs  in  a  Mo  +  1%  Ni  specimen  at  1400  °C,  which  lies  in  the  single-phase  region  where  the  bulk  liquid 
is  not  yet  stable,  (c)  Computed  thickness  of  stable  IGFs  vs.  the  Ni  bulk  fraction  at  1400  °C,  where  the  “X”  indicates 
the  Mo  +  1%  Ni  specimen.  Reprinted  from  an  AFOSR  supported  publication  in  Appl.  Phys.  Lett.  [9]. 
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Ni-doped  Mo:  Summary 

Shi  &  Luo,  APL  94:  251908  (2009) 


Fig.  13:  SEM  images  reveal  that  (a)  the  crystalline  5-NiMo  phase  does  not  wet  Mo  GBs,  but  (b)  the  Ni-rich  liquid 
phase  does  wet  Mo  GBs  completely,  (c)  HRTEM  and  AES  reveal  the  thermodynamic  stabilization  of  ~l-nm-thick 
quasi-liquid  IGFs  in  the  single-phase  region  (where  the  bulk  liquid  is  no  longer  stable).  Images  reprinted  from  an 

AFOSR  supported  publication  in  Appl.  Phys.  Lett.  [9]. 


Fig.  14:  Snapshots  of  the  atomistic  simulation  illustrate  that  (b)  a  GB  in  Ni-doped  Mo  is  significantly  more 
disordered  and  wider  than  that  in  (a)  pure  Mo.  (c)  Simulation  shows  that  GB  diffusivity  is  significantly  increased 
with  Ni  doping.  Figures  are  courtesy  of  our  collaborator,  Professor  Hao  Zhang,  Univ.  Alberta  in  Canada. 


Furthermore,  we  constructed  a  GB  diagram  for  Ni-doped  Mo,  which  is  shown  in  the  upper- 
right  corner  in  Fig.  15.  To  verify  the  correctness  of  this  computed  GB  diagram,  we  conducted 
well-controlled  sintering  experiments  to  extract  GB  diffusivity  as  a  function  of  the  temperature 
and  overall  composition.  The  experimental  procedure  and  detailed  results  are  documented  in 
Appendix  C.  These  experimental  results  can  be  well  explained  from  the  computed  GB  diagram. 
The  key  points  are  summarized  in  Fig.  15  and  discussed  below: 

•  Major  jumps  in  measured  GB  diffusivities  (labeled  by  the  red  stars  in  Fig.  15(A)  and  (B)) 
occur  at  a  “GB  solidus  line”  that  corresponds  to  a  computed  Tl  value  of  ~0.5  nm.  This  GB 
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solidus  line  (i.e.,  the  thick  red  dashed  line  in  Fig.  15),  where  presumably  quasi-liquid  IGFs 
start  to  fonn,  is  distinct  from  the  bulk  solidus  line.  Since  the  densification  represents  an 
average  effect  of  many  GBs,  a  gradual  transition  in  GB  diffusivity  was  expected  and 
observed. 

•  At  the  bulk  solidus  line  where  the  bulk  liquid  appears,  the  measured  (effective)  GB 
diffusivities  only  jump  moderately  (labeled  by  the  blue  stars  in  Fig.  15(A)  and  (E)). 

•  Below  the  bulk  peritectic  temperature,  there  is  no  discontinuity  in  the  measured  GB 
diffusivity  at  the  bulk  solvus  line  (labeled  by  the  light  purple  star  in  Fig.  15(B)). 

•  In  pure  Mo,  the  measured  GB  diffusivity  increases  with  increasing  temperature  (Fig. 

15(C)),  which  is  well  expected. 

•  In  the  single-phase  region,  the  measured  GB  diffusivity  appears  to  decrease  with  increasing 
temperature  at  a  constant  Ni  fraction  of  0.5  %  (Fig.  15(D)).  This  counterintuitive  result  can 
be  fully  explained  from  the  computed  GB  diagram.  This  abnormal  behavior  of  decreasing 
GB  diffusivity  with  increasing  temperature  is  due  to  the  fact  that  Mo-Ni  has  a  “retrograde” 
phase  diagram.  With  increasing  temperature,  this  fixed  composition  of  0.5  %  moves  away 
to  the  metastable  extension  to  the  solidus  line  (i.e.,  the  black  dotted  line  in  the  computed 
GB  diagram  in  Fig.  15).  Consequently,  the  free  energy  penalty  for  forming  undercooled 
liquid  increases  with  increasing  temperature.  Consistently,  the  computed  AL  decreases  with 
increasing  temperature,  which  is  shown  in  the  inset  of  Fig.  15(D).  Thus,  a  GB  can 
“solidify”  with  increasing  temperature! 

The  above  comparisons  clearly  support  the  correctness  and  usefulness  of  our  thermodynamic 
models  and  the  computed  GB  diagram. 


Bulk  Diffusion 


_ ►  T(°  C)  _ ►  r(°C)  _ _  7*(°C) 

Fig.  15:  Measured  GB  diffusivities  (panels  (A)  to  (E))  are  compared  to  the  computed  GB  diagram  for  Ni-doped  Mo 
(shown  in  the  upper-right  corner)  to  validate  the  correctness  and  usefulness  of  our  thermodynamic  models. 
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3.5.  Complexion:  A  Transformative  Concept  &  Contributions  via  This  Program 

Tang  et  al.  [32,  35]  suggested  an  extension  of  the  Gibbs  definition  of  bulk  phases  to 
equilibrium  GB  features  and  designated  them  as  "GB  complexions".  They  predicted  first-order 
and  continuous  GB  transitions  [32,  35],  resulting  in  abrupt  changes  in  GB-sensitive  transport, 
mechanical  and  physical  properties. 

In  2007,  Dillon  and  Harmer  (Lehigh  University)  observed  six  distinct  GB  complexions  in 
ALO3  (Fig.  16a)  [17,  36].  These  six  GB  complexions  can  be  grouped  into  three  classes: 

•  Complexion  I  (monolayer  adsorption),  Complexion  II  (clean/crystalline  GBs)  and 
Complexion  VI  (complete  wetting)  are  well  known. 

•  Complexion  V  (1-2  nm  IGFs),  discussed  in  §2.3,  are  also  widely  observed. 

•  Complexion  III  (bilayers)  and  Complexion  IV  (trilayers)  can  be  considered  to  be 
derivative  complexions  of  nanoscale  IGFs  (in  the  Pi’s  opinion). 

Dillon  et  al.  used  the  concept  of  GB  complexions  to  solve  an  outstanding  scientific  problem 
regarding  the  origin  of  abnonnal  grain  growth  [17]  and  "demystify  the  role  of  sintering  aids" 
[16].  They  further  proposed  a  provocative  concept  of  the  “kinetic  engineering  of  materials  via 
utilizing  GB  complexions”  [17]. 

In  our  AFOSR-YIP,  we  proposed  a  phenomenological  model  [3,  4,  37]  for  Dillon-Harmer 
complexions  from  an  analogy  to  the  colloidal  theory  (Fig.  16b).  The  basic  concepts  can  be 
explained  as  follows.  By  reducing  the  temperature,  or  dopant  activity,  a  quasi-liquid  IGF  should 
become  thinner  as  opposing  grains  experience  attractive  forces.  As  the  film  thickness 
approaches  atomic  sizes,  an  oscillation  in  excess  free  energy  versus  film  thickness  arises, 
resulting  in  free  energy  minima  when  the  thickness  corresponds  to  an  integer  of  atomic  sizes. 
Consequently,  additional  stable  GB  complexions  with  discrete  characteristic  thicknesses  form. 


o  Oxygen  Anion 


•  Dopant  Cation  •  Aluminum  Cation 


(h  ->  -foo) 
bulk  liquid 

VI 


GB  Solidus 


Bulk  Solidus 


(a)  (b) 

Fig.  16:  (a)  Schematic  illustrations  of  six  Dillon-Harmer  Complexions  in  A1203  from  Ref.  [17]: 
Complexion  I  -  the  (Langmuir-McLean)  monolayer  adsorbates;  Complexion  II  -  clean/crystalline  GBs;  Complexion 
III  -  adsorption  bilayers;  Complexion  IV-  adsorption  trilayers;  Complexion  V —  1-2  nm  thick  IGFs;  and 

Complexion  VI  -  completely  wetting  fdms. 

(b)  A  thermodynamic  model  developed  in  our  AFOSR-YIP  explaining  the  origins  of  the  Dillon-Harmer 
complexions  from  atomic  size  effects  [3,  4,  37].  Panel  (b)  is  adapted  from  an  AFOSR  supported  publication  in  Appl. 

Phys.  Lett.  [4]. 
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4.  Relevance  to  Aerospace  Materials  and  Broader  Impacts 


Comprehending  liquid-like  GB  “phases”  is  not  only  scientifically  challenging,  but  also 
technologically  important.  In  addition  to  critically  affecting  high-T  properties  and  the  materials 
fabrication  processing,  impurity-based  GB  “phases”  can  often  retain  upon  cooling  and  affect 
low-T  mechanical  and  physical  properties.  Table  2  lists  selected  technological  relevance. 


Table  2:  Selected  technological  roles  of  impurity-based  IGFs  or  analogous  liquid-like  GB  complexions.  This  table 
was  adapted  from  a  recent  critical  review  article  written  by  the  PI  [1], 


High-Temperature 
Material  Properties 

High-T  strength,  fracture  behaviors,  creep  resistance,  and 
oxidation  of  Si3N4-  and  SiC-  based  structural  ceramics 

[54-62] 

Superplasticity  of  Zr02-based  ceramics 

[63] 

WC-Co  composites 

[64] 

Refractory  metal  alloys 

[2,  30] 

Mechanical 
Properties  and 
Chemical  Stability 

Toughness,  strength,  fatigue,  and  wear  resistance  of  Si3N4,  SiC, 
and  A1203  based  structural  ceramics 

[65-71] 

Liquid  metal  embrittlement  for  Ni-Bi,  Cu-Bi,  and  Al-Ga 

[72-75] 

Corrosion  of  synroc 

[76] 

GB  embrittlement  of  W-Ni,  Cu-Bi,  and  other  systems 

[42] 

Solid-State  Activated 
Sintering 

ZnO-BL03;  Ce02-CoO 

[3,  4,  9,  44,  77,  78] 

W-Ni;  W-Pd;  W-Fe;  W-Co;  Mo-Ni 

[2,  30,  79] 

Ice-impurities 

[19,  20] 

Abnormal  Grain 
Growth 

Al203-(Y203+Si02);  Al203-(other  dopants) 

[17,31,36] 

Al-Ga,  W-Ni,  and  ZnO-BL03 

[80-83] 

Electronic  and 
Physical  Properties 

Ruthenate-based  thick-film  resistors 

[84] 

ZnO-BL03  based  varistors 

[29,  44,  85] 

(Ba,  Sr)Ti03-based  sensors  and  actuators 

[86-91] 

Thermal  conductivity  in  AIN 

[92,  93] 

High-7].  superconductors 

[94] 

Controlling  Microstructure:  GB  Diffusion  vs.  Migration: 

Nanoscale  IGFs  often  control  densification  [3,  4,  9,  44,  77,  78].  The  presence  of  impurity- 
based  liquid-like  GB  structures  can  also  promote  GB  migration.  This  mechanism,  which  differs 
from  the  classical  solute-drag  model,  has  been  observed  for  doped  ABCfi  [17,  31,  36],  ZnO- 
BBO3  [83],  W-Ni  [82],  and  Al-Ga  [80,  81].  Dillon  and  Hanner  [17,  36]  further  suggested  that 
the  existence  of  two  or  more  GB  complexions  with  significantly  different  GB  mobilities  will  lead 
to  detrimental  abnormal  grain  growth.  Practically,  it  is  often  desirable  to  achieve  full 
densification  with  minimum  grain  growth  (to  make  nanocrystalline  or  fine-grained  materials)  via 
pressure-less  sintering  (instead  of  using  "brute  force  methods").  Here,  we  recommend  using  GB 
diagrams  to  predict  a  kinetic  window  where  GB  diffusion  is  significant,  but  GB  migration  is 
limited.  This  supports  our  long-range  scientific  goal  of  developing  “GB  diagrams”  as  a  new 
tool  for  mechanism-informed  materials  design. 

High-Temperature  Structural  Materials: 

Quasi-liquid  IGFs  appear  more  frequently  in  high-temperature  structural  ceramics  and 
refractory  metals,  because  their  fonnation  will  reduce  the  otherwise  immense  GB  energies.  In 
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general,  the  presence  of  impurity-based  IGFs  reduces  the  strength  and  creep  resistance,  but 
increase  fracture  toughness.  Specifically: 

•  It  is  now  well  known  that  the  control  of  IGFs  is  critical  for  achieving  low-7"  toughness  and 
yet  optimizing  high-temperature  creep  resistance  for  Si3N4  based  ceramics  [65,  67,  68], 

•  For  SiC  based  ceramics,  a  recent  study  has  demonstrated  that  drying/crystallizing  IGFs  is 
possible,  leading  to  increased  creep  resistance  at  the  expense  of  fracture  toughness  [95]. 

•  IGFs  also  play  critical  roles  in  the  mechanical  properties  and  erosive  wear  behaviors  of 
AFO3  [70,  71]  and  the  superplasticity  of  Y2O3  stabilized  ZrC>2  [63]. 

•  While  the  presence  of  IGFs  generally  reduce  creep  resistance,  the  fonnation  of  Complexion 
I  in  (at  least)  Y2C>3-doped  AI2O3  significantly  increases  creep  resistance  [96]. 

•  IGFs  exist  in  WC-Co  composites  and  affect  their  sintering  and  mechanical  properties  [64]. 

•  IGFs  with  similar  characteristics  exist  at  hetero-phase  interfaces  (e.g.,  AkC^-TiN  [97],  BN- 
SiC  [98],  AFC^-metals  [99]).  They  affect  the  mechanical  properties  of  these  composites. 

Analogous  IGFs  in  binary  refractory  metals  (studied  in  this  AFOSR-YIP)  and  in  Mo-Si-B  based 
multicomponent  alloys  (to  be  investigated  in  our  new  AFOSR  project)  serve  as  relatively  simple 
model  systems  that  can  be  used  to  establish  and  validate  the  basic  high-temperature  interfacial 
thermodynamic  model.  In  future  studies,  we  plan  to  extend  the  model  to  more  complex  ceramic 
materials  (where  electrostatic  and  vdW  forces  may  complicate  interfacial  behaviors).  We  have  a 
particular  interest  in  analyzing  ultra  high  temperature  ceramics  (UHTC). 

GB-Controlled  Oxidation,  Corrosion  and  Embrittlement: 

The  presence  of  quasi-liquid  IGFs  may  critically  impact  the  GB-controlled  oxidation  and 
corrosion  resistance.  Recent  studies  have  revealed  the  important  role  of  nanoscale  IGFs  in  the 
liquid  metal  embrittlement  process  [72-75].  Furthennore,  impurity-based  IGFs  can  often  retain 
upon  cooling.  For  metallic  systems,  they  typically  result  in  GB  embrittlement.  Here,  using 
(abrupt)  first-order  GB  transitions  to  design  heat  treatment  procedures  to  alleviate  such 
detrimental  effects  can  be  a  potentially  transformative  approach. 

Functional  Ceramics: 

IGFs  can  critically  affect  the  tunable  conductivities  for  ruthenate  based  thick-film  resistors 
[84],  non-linear  I-V  characteristics  for  ZnO-BFCh  based  varistors  [29,  85],  functions  of  (Sr, 
Ba)TiC>3  based  sensors  and  actuators  [86-91],  thermal  conductivity  of  AIN  substrates  [92,  93], 
and  critical  currents  of  high  Tc  superconductors  [94].  These  functional  ceramics  also  have  a 
wide  range  of  applications  in  the  aerospace  industry. 
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5.  Executive  Summary  (Conclusions): 

Impurity-based  intergranular  films  or  IGFs  have  been  widely  observed  in  structural  and 
electronic  ceramics,  where  they  often  control  sintering,  microstructural  evolution  (nonnal  and 
abnonnal  grain  growth),  creep  resistance,  strength,  fracture  toughness,  and  some  electronic 
properties.  Under  the  support  of  this  three-year  AFOSR  Young  Investigator  Program,  we 
observed  and  systematically  studied  analogous  quasi-liquid  IGFs  in  W  and  Mo  based  binary 
refractory  alloys. 

Combined  experimental  and  modeling  studies  of  both  W  and  Mo  based  systems 
unequivocally  solved  an  outstanding  scientific  problem  regarding  the  origin  the  solid-state 
activated  sintering  in  refractory  metals,  which  had  puzzled  the  materials  community  for  more 
than  half  a  century.  This  study  showed  that  the  subsolidus  enhanced  densification  is  due  to  the 
increased  mass  transport  in  impurity-based  liquid-like  IGFs.  These  quasi-liquid  IGFs  are 
thennodynamic  stabilized  below  the  bulk  solidus  line,  and  their  formation  can  be  greatly 
promoted  via  minor  addition  of  certain  dopants  (sintering  aids).  For  the  first  time,  quantitative 
predictions  for  the  effectiveness  of  solid-state  sintering  aids  have  been  made  via  a  high- 
temperature  interfacial  thermodynamic  model.  This  model  quantitatively  predicted  the  onset 
subsolidus  activated  sintering  temperatures  for  Pd,  Ni,  Co,  Fe,  and  Cu  doped  W.  It 
systematically  explained  all  major  observations  about  activated  sintering  in  Ni-doped  Mo. 

CalPhaD  (Calculation  of  Phase  Diagrams)  methods  were  combined  with  a  statistical 
thennodynamic  model  (a  Miedema-type  “macroscopic  atom”  model)  to  predict  the 
thennodynamic  stability  of  nanometer-thick  undercooled  liquid  films  sandwiched  at  grain 
boundaries  (GBs).  This  model  predicted  that  a  “liquid  precursor”  could  appear  at  GBs  at  as  low 
as  60-85%  of  the  bulk  solidus  lines  in  certain  systems  (such  as  Pd  or  Ni  doped  W).  This  model 
was  further  refined  and  validated  via  comparing  with  various  experimental  data.  Furthermore,  a 
phenomenological  model  was  developed  to  consider  atomic  size  effects,  and  this  model 
produced  a  series  of  discrete  “GB  phases”  with  character  similar  to  those  six  “complexions” 
observed  by  Dillon  and  Harmer. 

We  further  proposed  a  long-range  scientific  goal  of  developing  “GB  diagrams”  as  a  new 
materials  science  tool.  Such  GB  diagrams  can  be  used  for  1)  optimizing  fabrication  pathways  to 
utilize  desired  GB  “phases”  during  processing  to  control  microstructural  evolution;  2)  designing 
heat  treatment  recipes  to  adjust  GB  structures  to  improve  mechanical  or  physical  properties;  and 
3)  forecasting  GB-related  high-temperature  properties.  Alternatively,  failing  to  appreciate  the 
presence  of  liquid-like  GB  structures  at  high  temperatures  can  result  in  a  significant  variability  in 
material  fabrication  and  properties  that  may  disqualify  their  successful  aerospace  applications. 

Ni-doped  Mo  was  selected  for  systematical  evaluation  via  characterizing  well-quenched 
specimens  and  thennodynamic  modeling.  In  contrast  to  prior  reports,  we  showed  that  the  5- 
NiMo  phase  does  not  wet  Mo  GBs  in  the  solid  state.  In  the  solid-liquid  two-phase  region,  the 
Ni-rich  liquid  wets  Mo  GBs  completely.  Furthermore,  atomic-resolution  HRTEM  and  AES 
characterization  of  well-quenched  specimens  unequivocally  demonstrates  that  nanometer-thick 
quasi-liquid  films  persist  at  GBs  into  the  single-phase  region  where  the  bulk  liquid  phase  is  no 
longer  stable.  To  verify  the  correctness  of  a  computed  GB  diagram  for  Ni-doped  Mo,  well- 
controlled  sintering  experiments  were  conducted  to  extract  GB  diffusivity  as  a  function  of  the 
temperature  and  overall  composition.  These  experimental  results  can  be  well  explained  from  the 
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computed  GB  diagram.  Specifically,  major  jumps  in  GB  diffusivities  occur  at  a  distinct  “GB 
solidus  line”.  In  contrast,  there  are  only  moderate  jumps  in  the  effective  GB  diffusivities  at  the 
bulk  solidus  line  where  the  bulk  liquid  appears,  and  there  is  no  jump  in  measured  GB  diffusivity 
at  the  bulk  solvus  line  below  the  bulk  peritectic  temperature.  In  pure  Mo,  the  measured  GB 
diffusivity  increases  with  increasing  temperature,  which  is  well  expected.  However,  with  the 
addition  of  0.5%  Ni,  the  measured  GB  diffusivity  appears  to  decrease  with  increasing 
temperature  (although  this  composition  lies  in  the  single-phase  region  for  all  temperatures).  This 
counterintuitive  result  can,  in  fact,  be  fully  explained  from  the  computed  GB  diagram.  The 
modeling  and  experimental  results  illustrated  that  a  GB  can  “solidify”  with  the  increasing 
temperature  at  certain  condition! 

The  formation  of  liquid-like  GB  structures  is  particularly  important  for  high-temperature 
materials  such  as  refractory  metals  and  ceramics.  This  is  because  the  GB  energy  typically  scales 
with  the  melting  temperature;  greater  GB  energies  can  stabilize  quasi-liquid  IGFs  to  larger 
undercooling  regions,  thereby  have  greater  impacts  on  microstructural  evolution  and  high- 
temperature  properties. 

The  feasibility  and  necessity  of  developing  GB  diagrams  as  new  materials  science  tool  for 
controlling  microstructural  evolution  and  forecasting  high-temperature  properties  are  proven  by 
this  AFOSR  program  on  binary  refractory  alloys.  There  is  a  critical  need  to  extend  the  developed 
high-temperature  interfacial  thennodynamic  models,  as  well  as  the  methods  to  compute  GB 
diagrams,  to  multicomponent  alloys  to  be  practically  useful.  These  models  and  methods  should 
be  further  extended  to  more  complex  ceramic  materials  in  future  studies. 


A  List  of  People  Involved  in  the  Research  Effort 

•  Jian  Luo,  the  PI  and  project  director  (who  received  1  summer  month  per  year  support  and 
worked  on  modeling  and  design  of  experiments) 

•  Xiaomeng  (Jimmy)  Shi,  a  Ph.D.  student  (who  received  full  support  and  worked  solely  on 
this  AFOSR  project  on  modeling  and  Mo-Ni  experiments) 

•  Jingou  Yin,  a  Ph.D.  student  (who  received  partial  support  and  worked  on  sintering  of  W) 

•  Archana  Kayyar,  a  Ph.D.  student  (who  received  partial  support  and  worked  on  a  regular 
solution  model  for  GBs  in  binary  alloys) 
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A  Special  Highlight  (Student  Award): 

Jimmy  Shi,  a  Ph.D.  student  supported  by  this  AFOSR  grant,  was  one  of  three  students 
selected  worldwide  for  a  Diamond  Award  for  2009  Graduate  Excellence  in  Material  Sciences 
(GEMS)  awards  from  the  American  Ceramic  Society,  Basic  Science  Division. 

The  GEMS  awards  are  selected  following  a  two-step  process  (written  applications  and  oral 
presentations  for  the  selected  finalists  at  the  annual  Materials  Science  and  Technology 
conference).  Diamond  awards  are  the  highest  catalog  for  the  GEMS  awards. 

Please  see  the  official  awards  page:  http ://ceramics .org/ acers-blog/2009-gems-awards/  as 
well  as  the  following  announcement  in  the  2010  January/February  issue  of  The  American 
Ceramic  Society  Bulletin. 


Diamond  Awards 


Jimmy  Shi 

Segregation.  Disordering 
and  Wetting  at 
Molybdenum  Grain 
Boundaries  and  Their 
Implications  in  Activated 
Sintering 

Clemson  University 


The  Americon  Ceramic  Society  Bulletin 
89(1).  pp.  9  (Jan/Feb  2010) 


2009  Graduate  Excellence  in  Materials 


Science  Awards  (Presented  by  ACerS  Basic  Science  Dlv.) 

Diamond  Awards 

Jimmy  Shi 


Segregation,  Discrderrtg 
and  Weang  at 
Molyhdenum  Gran 
Boundaries  and  Thar 
Implications  In  Activated  j 
Snoring 

Ctemson  University 

Hamid  Aatzl-Aitzemmi 
Phase  Field  Modelngol 
A  us  tent  e  Formation  Iron 

_  Ferine  CatUa 

Afaonni  Aggregates 

University  ot  British 
Columbia 

Craig  Przybyta 

Microstructure-  Sens  tire 
Extreme  Value 
Probatxlties  at  Fatigue  In 
Prrybyia  Adranced  Engineering 
Aloys 

Georgia  Institute  of 
Technology 


Sapphire  Awards 

HSudhakar  She) 

Bandgap-Reduced 
P-  Type  ZnO  Ft  ms  try 
Co-doping  at  Cu  and  6a 
ShM  tor  Imprortig 

Pholoelectro-  me  me  at 
Response 


Sapphire  Awards 

Weronika  Walkou 

menaces  ot  Site  on 
tttrtde  Ceramcs.  Direct 
Imgtng  and  Firs  l- 
Woi»„  Pmcptes  Calculations 
University  ot  Knots  at 
Chicago 


Yaodong  Yang 

Sold  Stale  Synthesis  ot 

Perorskjte-Spmet 

Nanocompostes 


Yang  Virginia  Tech 


C*l*ard 


Shew  Gittard 
Two-Photon 
Polymeric* t» i  ot 
Organically  Mod! led 
Ceramic  Materials 
North  Carolina  State 
University 

Sarthak  Haveila 

Noeel  Approaches  to  tie 
Growth  ot  Eptaanl 
BiMnO,  Thn  Pints 
Carnegie  Melon  University 


New  Jersey  institute  ot 
Technology 
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A  list  of  publications  that  acknowledged  this  AFOSR  grant: 

•  X.  Shi  and  J.  Luo,  "Grain  Boundary  Wetting  and  Prewetting  in  Ni-Doped  Mo,"  Applied 
Physics  Letters,  94,  251908  (2009). 

•  J.  Luo  and  X.  Shi,  "Grain  Boundary  Disordering  in  Binary  Alloys,"  Applied  Physics  Letters, 
92,  101901  (2008). 

•  J.  Luo,  "Subsolidus  Activated  Sintering  and  Liquid-like  Interface  Complexions"  (invited 
review  article),  Current  Opinion  in  Solid  State  and  Materials  Science,  12,  81-88  (2008) 

•  J.  Luo,  "Stabilization  of  Nanoscale  Quasi-Liquid  Interfacial  Films  in  Inorganic  Materials: 

A  Review  and  Critical  Assessment"  (invited  review  article,  image  featured  in  the 
journal  cover),  Critical  Reviews  in  Solid  State  and  Materials  Sciences,  32,  67-109  (2007) 

•  V.  K.  Gupta,  D.  H.  Yoon,  H.  M.  Meyer  III,  and  J.  Luo,  "Thin  Intergranular  Films  and 
Solid-State  Activated  Sintering  in  Nickel-Doped  Tungsten,"  Acta  Materialia,  55,3131- 
3142  (2007) 

•  J.  Luo,  S.  J.  Dillon  and  M.P.  Harmer,  “Interface-Stabilized  Nanoscale  Quasi-Liquid  Films” 
(invited  overview),  Microscopy  Today,  17  [4],  22-26,  July  2009. 

•  Z.  Wei,  G.  Zhang,  H.  Chen,  J.  Luo,  R.  Liu,  and  S.  Guo,  "A  Simple  Method  for  Evaluating 
Elastic  Modulus  of  Thin  Films  by  Nanoindentation,"  Journal  of  Materials  Research,  24 
[3],777-791  (2009). 

•  R.  H.  French,  V.  A.  Parsegian,  R.  Rodgomik,  R.  F.  Rajiter,  A.  Jagota,  J.  Luo,  D.  Asthagiri, 
M.  K.  Chaudhury,  Y.-M.  Chiang,  S.  Granick,  S.  Kalinin,  M.  Kardar,  R.  Kjellander,  D.  C. 
Langreth,  J.  Lewis,  S.  Lustig,  D.  Wesolowski,  J.  Wettlaufer,  W.-Y.  Ching,  M.  Finnis,  F. 
Houlihan,  O.  A.  von  Lilienfeld,  C.  J.  van  Oss,  and  T.  Zemb,  "Long  Range  Interactions  in 
Nanoscale  Science"  (invited  review  article),  Reviews  of  Modern  Physics,  in  press 

•  J.  Luo,  H.  Qian,  and  X.  Shi  "Interface  Stabilized  Nanometer-Thick  Disordered  Films,"  a 
review  chapter  in  the  2nd  edition  of  Encyclopedia  of  Nanoscience  and  Nanotechnology, 
edited  by  H.  S.  Nalwa,  American  Scientific  Publishers,  in  press. 

•  Xiaomeng  Shi,  “Grain  Boundary  Premelting  and  Activated  Sintering  in  Binary  Refractory 
Alloys,”  Ph.D.  dissertation,  Materials  Science  and  Engineering,  Clemson  University, 
Clemson,  South  Carolina,  USA,  May  2010 


A  list  of  presentations  that  acknowledged  this  AFOSR  grant: 

•  J.  Luo,  “Nanoscale  Quasi-Liquid  Interfacial  Films:  The  Interplay  of  Premelting,  Prewetting 
and  Multilayer  Adsorption,”  SYMPOSIUM:  Computational  Thennodynamics  and  Kinetics 
(organizers:  Jeffrey  J.  Hoyt,  Dallas  Trinkle),  2010  TMS  Annual  Meeting  &  Exhibition, 
Seattle,  WA,  USA,  February  14-18,  2010.  (Invited  talk) 

•  J.  Luo,  “Impurity-Based  Nanoscale  Quasi-Liquid  Interfacial  Films,”  Focused  Session  3: 
Computational  Design,  Modeling,  Simulation  and  Characterization  of  Ceramics  and 
Composites,  ICACC  10,  The  34th  International  Conference  and  Exposition  on  Advanced 
Ceramics  and  Composites,  Daytona  Beach,  Florida,  January  24-29,  2010.  (Invited  talk) 

•  J.  Luo,  “Thermodynamic  Stabilization  of  Subsolidus  Quasi-Liquid  Intergranular  Films,” 
SYMPOSIUM:  Structural  Transitions  and  Local  Defonnation  Processes  at  and  Near  Grain 
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Boundaries,  MS&T’09,  Materials  Science  &  Technology  2009  Conference,  Pittsburgh, 
October  25-29.  (Invited  talk) 

•  X.  Shi,  J.  Yin  and  J.  Luo,  “Segregation,  Disordering  and  Wetting  at  Molybdenum  Grain 
Boundaries  and  Their  Implications  in  Activated  Sintering,”  SYMPOSIUM:  Powder 
Metallurgy  Processing  and  Products,  MS&T’09,  Materials  Science  &  Technology  2009 
Conference,  Pittsburgh,  October  25-29.  (Student  presentation) 

•  J.  Luo,  “Liquid-Like  Interfacial  Complexions:  Grain  Boundaries  vs.  Free  Surfaces  and 
Metals  vs.  Ceramics,”  International  Workshop  on  Interfaces,  Bear  Creek,  PA,  October  20- 
23,  2009.  (Invited  talk) 

•  M.  P.  Hanner  and  J.  Luo,  “Grain  Boundary  Complexions,”  ICCCI-3,  The  Third 
International  Conference  on  the  Characterization  and  Control  of  Interfaces  for  High  Quality 
Advanced  Materials,  and  Joining  Technology  for  New  Metallic  Glasses  and  Inorganic 
Materials,  Kurashiki,  September  6-9,  2009.  (Invited  talk,  present  for  Prof.  Martin  Harmer) 

•  J.  Luo,  “Nanoscale  Liquid-Like  Interfacial  Films  and  Their  Roles  in  Controlling 
Microstructural  Evolution  and  High-Temperature  properties,”  poster  presentation,  AFOSR 
Workshop  on  Aerospace  Materials  for  Extreme  Environments,  St.  Louis  (Clayton), 
Missouri,  August  3-5,  2009. 

•  J.  Luo  and  X.  Shi,  “Grain  boundary  diagrams:  from  activated  sintering  to  materials-by- 
design,”  The  Role  of  Grain  Boundary  in  Sintering,  Symposium  1 :  Design,  Modeling  and 
Simulation  of  Ceramic  Interfaces,  The  8th  Pacific  Rim  Conference  on  Ceramic  and  Glass 
Technology  The  8th  Pacific  Rim  Conference  on  Ceramic  and  Glass  Technology, 
Vancouver,  June  2,  2009.  (Invited  talk) 

•  J.  Luo,  H.  Qian,  and  X.  Shi,  "Stabilization  of  Quasi-Liquid  Interfacial  Films:  Adsorption, 
Disordering  and  Wetting,"  Ceramic  Surfaces,  Grain  Boundaries  and  Interfaces,  Materials 
and  Technology  2008:  Fundamentals  and  Characterization,  Pittsburg,  PA,  Oct  5-9,  2008. 

(Invited  talk) 

•  J.  Luo,  "Interface-Stabilized  Nanoscale  Quasi-Liquid  Films,"  The  5th  International 
Conference  on  Advanced  Materials  and  Processing  (ICAMP-5),  Harbin,  China,  September 
2-5.  (Invited  talk) 

•  J.  Luo,  "  Free-Surface  Counterparts  and  Metallic  Counterparts  to  Equilibrium  Intergranular 
Films  in  Ceramics:  Two  Related  Studies,"  Gordon  Research  Conference,  Solid  Studies  in 
Ceramics,  Proctor  Academy,  New  Hampshire,  August  10-15,  2008.  (Gordon  Conference 
Invited  Lecture) 

•  J.  Luo,  "Subsolidus  Quasi-Liquid  Interfacial  Films,"  International  Workshop  on 
Functionality  of  Interface  (jointed  sponsored  by  Max -Plank  Institution  and  Chinese 
Academy  of  Science  and  co-organized  by  Dr.  Hui  Gu  and  Dr.  Manfred  Ruhle),  April  16- 
17,  2008,  Shanghai,  China.  (Invited  talk) 

•  J.  Luo,  H.  Qian,  X.  Shi  and  V.  K.  Gupta,  "Stabilization  of  Nanometer-Thick  Subsolidus 
Quasi-Liquid  Interfacial  Films,"  oral  presentation,  Hume-Rothery  Symposium  -  Nanoscale 
Phases,  TMS  2008  Annual  Meeting  and  Exhibition,  New  Orleans,  LA,  March  10,  2008. 

(Invited  talk) 
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•  X.  Shi  and  J.  Luo,  "Segregation  and  Disordering  at  Tungsten  Grain  Boundaries  in  Well- 
Quenched  Specimens,"  poster  presentation,  Materials  and  Technology  2008:  Fundamentals 
and  Characterization,  Pittsburg,  PA,  Oct  5-9,  2008  (Student  presentation) 

•  X.  Shi  and  J.  Luo,  "Grain  Boundary  Wetting,  Prewetting  and  Activated  Sintering  in  Ni- 
Doped  Mo,"  poster  presentation,  Gordon  Research  Conference,  Solid  Studies  in  Ceramics, 
Proctor  Academy,  New  Hampshire,  August  10-15,  2008  (Student  presentation) 

•  J.  Luo,  "Nanoscale  Quasi-Liquid  Grain  Boundary  Films  in  Refractory  Metals,"  AFOSR 
High  Temperature  Aerospace  Materials  program  review  meeting.  May  12-16,  2008, 

Vienna,  VA. 

•  J.  Luo,  X.  Shi  and  V.  K.  Gupta,  "A  Quantitative  Model  for  Solid-State  Activated  Sintering 
and  Recent  Observations,"  oral  presentation,  The  32nd  International  Conference  & 
Exposition  on  Advanced  Ceramics  and  Composites,  Jan.  27-31,  2008 

•  X.  Shi  and  J.  Luo,  "Solid-State  Wetting  and  Activated  Sintering  in  Ni-Doped  Mo,"  oral 
presentation,  International  Symposium  on  Defects,  Transport  and  Related  Phenomena, 
Materials  and  Technology  2007:  Fundamentals  and  Characterization,  Detroit,  MI,  Sept  16- 
20,  2007.  (Student  presentation) 

•  X.  Shi  and  J.  Luo,  "A  Model  for  Grain  Boundary  Disordering  and  Solid-State  Activated 
Sintering  in  Doped  Tungsten,"  poster  presentation,  Gordon  Research  Conference,  Solid 
Studies  in  Ceramics,  Proctor  Academy,  New  Hampshire,  August  6-8  2007.  (Student 
presentation) 

•  J.  Luo,  "  Disordered  Interfacial  Films  and  Activated  Sintering:  Phenomenological 
Similarities  among  Ice,  Ceramics  and  Refractory  Metals,"  Refractory  Metals  2007,  2007 
TMS  Annual  Meeting  and  Exhibition,  Orlando,  FL,  February  28,  2007 

•  J.  Luo,  “Liquid-Like  Interfacial  Complexions,”  Karlsruhe  Institute  of  Technology,  Institute 
for  Ceramics  in  Mechanical  Engineering,  Karlsruhe,  Germany,  December  9,  2009.  (Invited 
seminar) 

•  J.  Luo,  “Impurity-Based  Quasi-liquid  Interfacial  Films:  The  Interplay  of  Premelting  and 
Prewetting,”  The  Materials  Science  Colloquium  Lecture,  the  Max  Planck  Institute  for 
Metals  Research,  Stuttgart,  Germany,  December  7,  2009.  (Invited  colloquium  lecture) 

•  J.  Luo,  “Equilibrium  Nanoscale  Interfacial  Films:  Disordering,  Adsorption  and  Wetting,” 
Departmental  Seminar,  Lehigh  University,  Lehigh,  PA,  February  24,  2009.  (Invited 
seminar) 

•  J.  Luo,  "Nanoscale  Quasi-Liquid  Interfacial  Films,"  Harbin  Institute  of  Technology, 

Harbin,  China,  September  5,  2008.  (Invited  seminar) 

•  J.  Luo,  "Interface-Stabilized  Nanometer-Thick  Disordered  Films,"  University  of  South 
Carolina  NanoCenter,  April  4,  2007.  (Invited  seminar) 

•  J.  Luo,  "Interface-Stabilized  Nanoscale  Quasi-Liquid  Films:  from  ice  skating  to  materials 
by  design,"  Mechanical  Engineering  Graduate  Seminar  Series,  Louisiana  State  University, 
Baton  Rouge,  March  23,  2007.  (Invited  seminar) 
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Appendices 


The  following  three  appendices  are  adapted,  in  part,  from  the  relevant  chapters  in  the 
doctoral  dissertation  of  Xiaomeng  (Jimmy)  Shi,  entitled  “Grain  Boundary  Premelting  and 
Activated  Sintering  in  Binary  Refractory  Alloys”  (Ph.D.,  Materials  Science  and 
Engineering,  May  2010,  Clemson  University,  Clemson,  South  Carolina,  USA). 
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Appendix  A: 

Experimental  Procedures 


(A.l)  Raw  Materials 

High  purity  molybdenum  (99.999%)  with  an  average  particle  size  of  20  pm  was  purchased 
from  Alfa  Aesar.  To  achieve  a  homogeneous  distribution  of  nickel  in  molybdenum,  nickel 
chloride  (NiCl2'6H20,  Alfa  Aesar)  was  mixed  with  molybdenum  in  a  solution/suspension.  Then, 
the  slurry  was  dried  in  an  oven  at  90  °C.  Afterwards,  the  powders  were  annealed  in  a  tube  furnace 
at  600  °C  for  1  hour  under  flowing  Ar  +  5%  H2  to  reduce  nickel  chloride  to  nickel. 

The  physical  properties  of  purchased  nickel  and  molybdenum  are  listed  in  Table  A.l.  The 
binary  phase  diagram  of  Ni-Mo  is  shown  in  Fig.  A.l.  We  investigated  the  GB  phenomena  and 
sintering  behaviors  in  three  phase  regimes:  the  single-phase  regime  (BCC),  the  sub-peritectic 
two-phase  regime  (BCC  +  8-NiMo),  and  the  solid-liquid  two-phase  regime  (BCC  +  liquid). 


Table  A.1  Physical  properties  of  nickel  and  molybdenum. 


Element 

Nickel  (Ni) 

Molybdenum  (Mo) 

Atomic  number 

28 

42 

Atomic  weight 

58.69 

95.94 

Density  (20  °C  g/cm3) 

8.90 

10.28 

Crystal  structure 

FCC 

BCC 

Atomic  radius  (nm) 

0.125 

0.139 

Ionic  radius  (nm) 

0.069 

0.062  (+6e),  0.070  (+4e) 

Most  common  valance 

2+ 

4+,  6+ 

Melting  point  (°C) 

1455 

2623 

Atomic  Percent  Molybdenum 


Weight  Percent  Molybdenum 


Figure  A.1  Phase  diagram  of  the  Ni-Mo  system.  The  inset  is  an  enlarged  view  of  the  regimes 

where  we  conducted  experiments. 
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(A.2)  Specimens  Preparation 


Ni-Mo  green  specimens  were  prepared  using  a  standard  powder  metallurgy  method.  The 
calcined  Ni-doped  Mo  powders  were  grinded  gently  using  a  spatula  or  mortar.  Thereafter, 
cylinder-shape  specimens  were  prepared  at  280  MPa  using  a  hand  press.  The  pressed  green 
specimens  had  a  diameter  of  ~6.36  mm  and  an  average  height  of  8-9  mm.  To  prevent  the 
formation  of  a  transient  liquid  and  to  ensure  the  chemical  homogeneity,  all  green  specimens  were 
pre -sintered  at  1100  °C  for  12  hours  in  flowing  Ar  +  5%  H2.  Then,  the  specimens  were  placed  in 
an  oil-quench  furnace  at  a  desired  temperature  and  isothermally  sintered  for  3  hours.  Afterwards, 
the  specimen  was  dropped  into  room-temperature  silicone  oil  in  less  than  1  second  inside  the 
quench  chamber  under  vacuum.  The  well-quenched  specimens  were  characterized  by  a  scanning 
electron  microscope  (SEM,  Hitachi  4800)  and  a  HRTEM  (Hitachi  9500,  300  kV). 

For  sintering  experiments,  the  pre -sintering  was  conducted  at  1000  °C  for  1  hour  and 
isothermal  sintering  was  carried  out  in  another  quench  furnace.  The  lower  pre-sintering 
temperature  and  shorter  pre -sintering  time  enable  us  to  have  a  bigger  window  to  extract  the  GB 
diffusivity  via  densification  data. 


(A.3)  Oil-Ouench  Furnace  and  Quenching  Experiments 

A  high-temperature  quench  furnace  was  built,  which  is  shown  in  Fig.  A.2.  Only  one 
specimen  may  be  placed  on  the  sample  holder  for  this  furnace  for  each  run.  The  sample  holder  is 
comprised  of  an  alumina  tube  and  a  spiral-shape  Mo  wire  container  (Fig.  A.3).  After  isothermal 
sintering  was  completed,  we  pulled  the  quench  handle  to  let  the  specimen  to  fall  into  quench  oil 
(diffusion  pump  silicone  oil),  which  should  happen  within  approximately  1  second. 


Figure  A.2  (a)  Schematic  illustration  and  (b)  picture  of  the  oil-quench  furnace. 
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(A.4)  Gas-Quench  Furnace  and  Sintering  Experiments 

We  assembled  a  second  furnace  for  specially  designed  sintering  experiments  to  extract  GB 
diffusivity  as  a  function  of  temperature  and  overall  composition  potential.  A  schematic  of  this 
furnace  is  shown  in  Fig.  A.4.  It  is  a  vertical  tube  furnace.  When  the  furnace  was  running,  Ar  +  5% 
H2  gas  is  flowing  through  the  system  to  protect  specimens  from  oxidation.  This  furnace  is  unique 
in  that  a  mechanism  of  loading/unloading  samples  was  built  so  that  the  specimens  can  be  inserted 
into  the  hot  zone  of  the  furnace  after  the  furnace  reaches  pre-set  temperature.  The  specimens  can 
also  be  taken  out  for  gas  quench  after  sintering  by  pulling  out  the  crucible  back  to  the  cold  zone. 
Thus,  the  effects  of  heating  and  cooling  ramps  are  minimized.  In  all  experiments,  a  second 
thermocouple  was  placed  in  the  crucible  to  monitor  the  actual  temperature. 

The  density  of  the  sample  was  determined  by  measuring  the  dimension  and  mass.  The  green 
specimens  were  cylindrical,  but  these  specimens  deformed  to  truncated  cone  shape  after  sintering. 
Therefore  three  measurements  were  conducted  at  different  heights  to  calculate  the  volume. 


Figure  A.4  Schematic  illustration  of  the  gas-quench  furnace. 
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(A.5)  Extraction  of  Diffusivity  Data 


Through  measuring  the  density  (and  grain  size)  as  a  function  of  time  during  sintering,  we  can 
calculate  the  densification  rates.  Then,  the  bulk  or  GB  diffusivity  can  then  be  determined  by 
fitting  the  well  established  sintering  models.  The  relevant  sintering  models  are  described  here. 
According  to  Coble,  densification  is  divided  into  three  stages:  the  initial  stage,  the  intermediate 
stage,  and  the  final  stage  (Fig.  A.5)  [7].  The  relevant  sintering  models  are  discussed  below. 


Figure  A.5  Typical  relative  density  vs.  time  during  sintering  [7]. 


Figure  A.6  Models  for  (a)  initial,  (b)  intermediate,  and  (c)  final  stage  sintering  [7]. 
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Initial-Stage  Sintering: 

The  initial  stage  sintering  is  often  modeled  using  a  two-sphere  model  as  shown  in  Fig.  A. 6(a). 
The  Johnson  model  [5]  is  the  most  widely  used  for  the  initial  stage  of  sintering.  The  key 
assumptions  are: 

•  Transport  due  to  dislocation  or  small  angle  grain  boundary  is  negligible; 

•  The  grain  geometry  can  be  represented  by  two  round  grains  in  contact  with  each  other; 

•  The  surface  and  GB  energy  are  isotropic; 

•  Vacancy  concentration  potentials  are  equal  for  bulk  diffusion  and  GB  diffusion;  and 

•  The  densification  rate  is  small  (AL/L0  <  0.05). 

In  the  Johnson  model,  the  densification  vs.  time  is  represented  by: 


a l=  K  rP-D 

L0  [  kTGp 


(A- 1 ) 


where: 


K  =  (A.l-6.9)*8;  m  =  -0.4-0. 5,  p  =  3 ,D  =  Dxl ;  for  bulk  diffusion  dominant  sintering; 

and 


K=  (1.6-3)*  16;  m  =  -0.3,  p  =  4 ,D  =  8 DGB  for  GB  diffusion  dominant  case. 

In  the  above  equations,  AL/L0  is  the  linear  shrinkage;  ys  is  the  surface  energy;  Q  is  the  volume  of 
single  vacancy;  Dxl  is  the  bulk  diffusivity;  DGB  is  the  GB  diffusivity;  8  is  the  GB  width  (which  is 
arbitrarily  selected  to  be  1  nm  in  most  prior  studies);  k  is  the  Boltzmann  constant;  T  is  the 
sintering  temperature;  and  G  is  the  grain  size.  German  also  proposed  a  similar  model  [3,6]. 


Intermediate-Stage  Sintering: 

In  the  intermediate  stage  (Fig.  A.6  (b)),  a  continuous  interconnected  pore  channel  forms  and 
shrinkage  is  large,  but  grain  growth  is  expected  to  be  mild  because  GBs  are  pinned  by  the  pore 
channels.  Coble  developed  a  model  for  intermediate-stage  sintering  [1],  with  the  following  key 
assumptions: 

•  GBs  are  the  sink  of  vacancies  and  pores  are  the  source  of  vacancies; 

•  The  grain  geometry  is  represented  by  a  tetrakaidecahedra  (Fig.  A.6(b))  with  a  continuous 
pore  channel  at  triple -grain  junctions;  and 

•  The  pore  channels  have  cylindrical. 

If  the  intermediate  stage  sintering  is  controlled  by  bulk  diffusion,  the  porosity  (Pv)  and  its 
derivative  with  respect  to  time  can  be  expressed  as: 


P  -P°  = 


-336  Dm- 


/P 


kTG 


(A.2) 


-  =  -33  6^^ 
dt  kTG 3 


(A.  3) 


On  the  other  hand,  if  the  intermediate  stage  sintering  is  controlled  by  GB  diffusion,  the  porosity 
and  its  derivative  with  respect  to  time  can  be  expressed  as: 
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(A.  4) 


pl  -p4 = _128i  SDgbYJ^  t 


3 

dP1 


=  -1281 


kTG 

SDgbY^ 


dt  kTG 4 

The  above  equations  assume  a  constant  grain  size.  Here  Pv  is  porosity  and  it  is  defined  as: 

Pv=(]-prelanJ  (A.5) 

where  Prelative  is  the  relative  density. 


Final-Stage  Sintering: 

After  the  interconnected  pore  channels  break  up  and  pores  become  isolated,  the  sintering 
enters  the  final  stage  (Fig.  A.6  (c)).  Coble  also  developed  a  final-stage  sintering  model  [1].  This 
model  is  not  discussed  here  because  we  did  not  use  final-stage  sintering  in  our  experiments  to 
extract  GB  diffusivities. 


(A.6)  TEM  Specimen  Preparation 

HRTEM  was  employed  to  characterize  the  structure  of  GB.  To  prepare  ~3  nun  self- 
supporting  TEM  specimens,  we 

( 1 )  used  Buehler  low-speed  diamond  saw  to  cut  a  sintered  sample  into  thin  pellets  which  had 
the  same  diameter  as  the  sintered  sample  (~5.8  mm)  and  was  ~  1  mm  in  thickness; 

(2)  polished  the  cut  pellets  on  both  sides  on  emery  papers  of  grits  320,  400,  600,  800,  and 
1000  (in  a  sequence)  to  reduce  the  thickness  to  about  0.2  mm; 

(3)  used  a  South  Bay  Technology  (SBT)  abrasive  disc  cutter  to  cut  -3  nun  diameter  discs; 

(4)  further  polished  the  disc  (0.2  nun  thick  and  -3  mm  in  diameter)  on  emery  paper  of  grit 
1000  and  then  diamond  lapping  film/diamond  paste  with  9  pm,  3  pm  and  1  pm  particle 
sizes  (in  a  sequence),  to  produce  a  final  foil  of  -50  pm  thick; 

(5)  used  a  Fischione  (Model  150D)  precision  dimple  grinder  to  dimple  grind  the  50  pm  thick 
foil  on  both  side  with  a  dimple  of  about  22  pm  in  depth  until  the  thinnest  section  in  the 
center  was  about  a  few  microns;'  and 

(6)  employed  a  Fischione  Model  1010  low  angle  argon  ion  mill  to  further  thin  the  center  of 
the  foil  until  a  perforation  appear  in  the  foil  and  good  electron  transparency  was  achieved. 

For  ion-milling  process,  key  adjustable  parameters  include  accelerating  voltage,  milling 
current,  milling  angle  and  milling  time.  Table  A.2  listed  a  few  optimized  milling  recipes  that 
produced  good  TEM  foils  for  Mo-Ni  specimens.  We  used  the  Hitachi-9500  300  kV  HRTEM  to 
image  our  specimen.  If  the  electron  transparency  was  not  good  enough,  additional  ion-milling  can 


For  hard  metals  such  as  Mo  or  W,  it  is  recommended  to  reduce  thickness  to  ~50  pm.  For  other  soft  metals 
or  ceramics  with  high  ion-milling  etch  rates,  -100  pm  is  generally  recommended 
1  Although  people  generally  leave  about  20  pm  thick  foil  in  the  center  of  the  specimen  and  use  a  ion-mill  to 
perforate  it,  in  our  case  of  hard  metals  with  low  ion-milling  etching  rates,  we  must  keep  the  final 
thickness  in  the  center  of  the  specimens  as  thin  as  possible  to  reduce  milling  time. 
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be  carried  out  to  further  thin  the  specimen.  Usually,  good  HRTEM  images  can  be  obtained  if  one 
can  find  a  large  field  of  transparent  region  from  the  hole  in  the  center  to  the  peripheral  edge. 


Table  A.2  Representtaive  ion-milling  recipes  for  Ni-Mo  specimens. 


Recipe  #1 


Accelerating 

Voltage 

Milling  Current 

Milling  Angle 

Stage 

Temperature 

Milling  Time 

5  kV 

5.5  mA 

12° 

Liquid  N2 

90  mins 

5  kV 

5.5  mA 

9° 

Liquid  N2 

30  mins 

5  kV 

5.5  mA 

5° 

Liquid  N2 

30  mins 

Recipe  #2 


Accelerating 

Voltage 

Milling  Current 

Milling  Angle 

Stage 

Temperature 

Milling  Time 

5  kV 

5.5  mA 

12° 

Liquid  N2 

120  mins 

5  kV 

5.5  mA 

9° 

Liquid  N2 

60  mins 

5  kV 

5.5  mA 

5° 

Liquid  N2 

30  mins 

(A. 7)  Auger  Electron  Spectroscopy  and  Specimen  Preparation 

Auger  Electron  Spectroscopy  (  AES)  was  used  to  identify  elemental  composition  of  surfaces 
by  measuring  the  energies  of  Auger  electrons.  Compositional  depth  profiling  was  performed  by 
using  an  independent  ion  beam  to  sputter  the  specimen  surface  while  collecting  the  Auger 
electron  spectrum.  The  Auger  experiments  were  conducted  at  Oak  Ridge  National  Laboratory 
(ONRL)  via  the  Eligh  Temperature  Materials  Laboratory  (EITML)  user  facility  program  using  a 
PEU  680  SAM  in  collaboration  with  Dr.  H.  M.  Meyer  111. 

To  prevent  any  contamination,  the  specimens  were  fractured  in  situ  in  an  ultra  high  vacuum 
(UHV)  chamber.  Analysis  of  GBs  is  possible  because  the  Ni-doped  Mo  sintered  specimens  are 
known  to  undergo  brittle  intergranular  fracture.  A  compositional  depth  profiling  was  conducted. 
The  sputtering  rate  was  calibrated  to  be  2  nm  per  minute  for  a  Si02  standard  at  the  same 
sputtering  conditions. 

Fig.  A.7  shows  schematically  a  specimen  fixture  for  in  situ  fracture  in  the  Auger  chamber. 
The  procedure  to  make  AES  specimens  is  described  as  follows: 

(1)  A  Buehler  low-speed  diamond  saw  was  used  to  cut  ~1  mm  thick  disc  from  sintered  sample. 

(2)  The  same  diamond  saw  was  employed  to  cut  rectangular  specimens  of  ~1  mm  x  1  mm  x 
3.5-5  mm. 

(3)  Two  notches  were  made  using  this  diamond  saw  on  midsection  of  specimen. 

(4)  After  going  through  a  rigorous  clean  procedure  (described  below),  the  notched  specimen 
was  fixed  into  stainless  tube  on  both  ends  using  SEM-grade  conductive  epoxy. 

(5)  The  epoxy-fixed  specimen  was  cured  in  a  90  °C  oven  for  6-8  hours.  After  cool  down,  the 
specimen  was  wrapped  with  foam  carefully  and  packed  in  a  plastic  box. 

To  minimize  the  carbon  contamination  in  AES  specimen  introduced  by  the  cutting  fluid,  a 
special  clean  procedure  was  implemented  to  all  the  notched  specimens  before  making  epoxy 
fixture.  The  specimens  were 

( 1 )  ultrasonically  cleaned  with  acetone  for  20  mins  to  remove  all  the  polar  organic 
contamination  and  then  let  dry  fully; 
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(2)  ultrasonically  cleaned  with  toluene  for  20  mins  to  remove  all  the  non-polar  organic 
contamination  and  then  let  dry  fully; 

(3)  ultrasonically  cleaned  with  hexane  for  20  mins  to  remove  all  the  non-polar  organic 
contamination  and  then  let  dry  fully; 

(4)  ultrasonically  cleaned  with  ethanol  for  20  mins; 

(5)  soaked  in  ethanol  overnight  24  hours  to  remove  all  the  polar  organic  and  hydrophilic 
contamination;  and 

(6)  dried  in  a  90  °C  oven  for  3  hours  to  evaporate  all  the  residue  organic  contaminant. 


Screw  to  tighten 
/  the  sample 

Sample  holder 
base 

The  hole  in  the 
center  of  sample 
holder's  dimension: 
ID=3 ,65mm 
h=8.37mm 


Figure  A.7  Schematic  illustration  of  an  AES  specimen  fixture  for  in  situ  fracture. 


with  notch 


Tube  used  to 
support  W 
sample.  It  svlll  sit 
in  the  hole  of 
sample  bolder 
base 

For  the  tube: 
OD=-3.27nim 
ID=~1.2mni 
h=9-10mm 


(A. 8)  Grain  Boundary  Wetting  Experiments 

A  solid-state  wetting  experiment  was  performed  using  a  Mo  alloy  with  12.5  at.  %  Ni.  This  Ni 
content  represents  ~15  vol.  %  of  S-NiMo  compound  phase  at  the  sintering  temperature  (T=  1344 
and  1495  °C).  Ni  (99.999%;  Alfa)  and  Mo  (99.999%;  Alfa)  powders  were  mixed  in  distilled  de¬ 
ionized  water  and  dried.  The  mixture  was  compacted  in  a  die  at  280  MPa  to  produce  green 
specimen  with  ~  6.35  mm  in  diameter  and  8-9  mm  in  height.  Thereafter,  specimen  was 
presintered  at  1100  °C  for  12  hours  and  then  isothermally  sintered  at  a  desired  temperature  for  2 
hours  in  flowing  Ar  +  5%  H2.  The  sintered  specimens  were  polished  and  characterized  with  a 
Elitachi  S-4800  SEM  equipped  with  an  Oxford  EDX  detector. 


(A.9)  Grain  Size  Measurements 

Grain  size  evolution  of  specimen  was  analyzed.  A  Elitachi  S-4800  field  emission  SEM  and 
Olympus  optical  microscope  were  used  to  measure  the  grain  size.  The  specimens  were  polished 
to  mirror  finish  on  emery  papers  and  alumina  slurries  (or  diamond  paste).  Thereafter,  they  were 
etched  at  room-temperature  with  30%  volume  hydrogen  peroxide  to  make  grain  boundaries 
visible.  The  specific  procedure  for  producing  polished  and  etched  samples  for  grain  size  analysis 
is  described  as  follows: 

(1)  A  Buehler  low  speed  diamond  saw  was  used  to  cut  discs  from  sintered  pellets. 

(2)  The  specimens  were  then  roughly  polished  with  emery  papers  of  grit  sizes  of  240,  400,  600, 
800,  1000  in  a  sequence  (while  always  orienting  the  scratches  normal  to  the  direction  of 
prior  polishing). 
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(3)  Then,  specimens  were  finely  polished  with  alumina  suspensions  in  water  (or  diamond 
pastes)  to  get  a  mirror  finish. 

(4)  The  polished  Ni-Mo  specimens  were  etched  at  room-temperature  using  hydrogen  peroxide 
for  45  second  to  1  min. 
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Appendix  B: 
Thermodynamic  Models 


(B.l)  Phenomenological  Thermodynamic  Model 


A  phenomenological  thermodynamic  model  can  be  proposed  as  an  extension  to  a  premelting 
model  for  unary  systems  [1,2,3],  In  this  model  for  binary  alloys,  the  excess  GB  energy  (referred 
to  a  mixture  of  equilibrium  bulk  phases)  of  a  subsolidus  liquid-like  IGF  in  a  binary  A-B  alloy  is 
expressed  as: 

ax{h)  =  gs  -(fiArA+{iBrB)  =  2  rcl  + a  Gamorplfi  +  crinter6cial(/o ,  (b.i) 

where  Gs  is  excess  free  energy  according  to  the  Gibbs  definition,  //  is  the  bulk  potential,  /~is  the 
adsoiption  (GB  excess),  h  is  the  film  thickness,  yc\  is  the  interfacial  energy  of  the  crystal-liquid 
interface  (which  is  well  defined  when  the  two  interfaces  are  well  separated  or  h  +co),  AGamorph 
=  Giiquid  -  Gcrystai  is  the  volumetric  free  energy  for  forming  an  undercooled  liquid  (assuming  a 
uniform  film),  and  Ginterfaciai  is  the  interfacial  potential  (where  dom[crklcrfdh  is  the  well-known 
Derjaguin  disjoining  pressure).  The  interfacial  potential  (Ginterfaciai^))  represents  the  interactions 
of  two  interfaces  when  the  film  is  thin,  and  it  is  the  sum  of  all  short-  and  long-range  interfacial 
interactions  using  h  =  +co  as  the  reference  point.  By  definition, 

| ^nterfacial  (0)  =  (/gb  ~  2V c,  )  S  ~AX 

lfTnterfacial(+GO)  =  0  ’  (R2) 


where  we  define  yg°b  =  <7X  (0)  as  the  excess  free  energy  of  a  “dry”  GB,  which  is  in  general 
different  from  the  equilibrium  ygb .  The  equilibrium  GB  energy  ( ygb )  corresponds  to  the  global 
minimum  in  cr\h),  which  specifies  an  “equilibrium”  thickness: 


dcjfh) 

dh  h-K. 


(B3) 


One  may  further  define  a  dimensionless  interfacial  coefficient  as 

^"mterlacial  ^ ^  ) 

which  satisfies 


Ay 


ff(  0)  =  0 
l/(-H»)  =  1 


(B4) 


(B.5) 
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Then,  Eq.  (B.l)  can  be  rewritten  as: 

A  am  s  =  *r-m+*G.mp/t  (B 

A  liquid-like  IGF  of  thickness  h  can  be  thermodynamically  more  stable  than  a  dry  GB  if 
A  ex(/z)  <  0 .  As  such,  an  estimation  of  the  IGF  thickness  can  be  given  as: 


h 

eq  A  G 


■RK)< 


Ay 


amorph 


AG 


-X 


amorph 


(B.7) 


Here  the  subscript  “L”  is  used  to  denote  that  AGamorphis  calculated  using  a  reference  composition 
on  the  stable  or  metastable  liquidus  line.  (Possible  strategies  of  selecting  this  hypothetic 
reference  film  composition  are  discussed  in  §B.2,  among  which  this  assumption  leads  a 
theoretically  elegant  treatment  of  the  phenomenological  model.)  If  the  interfacial  coefficient 
takes  a  simple  exponentially  decaying  form  with  a  characteristic  coherent  length  of  c  (as  what  is 
commonly  assumed  for  the  premelting  theories  for  one-component  metals): 

f(h)«  \-eh'R  (B.8) 


minimizing  Acr(h )  leads  to: 


hEQ*{- w£/Z), 


(B.9) 


(B.2)  Through-Thickness  Gradients,  Film  Composition  and  Simplifications 

Through-thickness  compositional  and  structural  gradients  generally  exist  in  an  IGF. 
Correspondingly,  AGamorph  is  calculated  after  a  (somewhat  subjective)  selection  of  a  reference  film 
composition  (for  a  hypothesized  uniform  film  of  undercooled  liquid  as  a  reference  state);  the 
remaining  volumetric  free  energies,  along  with  the  excess  free  energies  associated  with  the 
compositional  and  structural  gradients,  is  then  be  considered  in  the  interfacial  potential  so  that  the 
thermodynamic  treatment  remains  rigorous.  Prior  studies  have  assumed  this  reference 
composition  to  be  the  average  film  composition  [4,5]  or  the  composition  that  maximizes  a 
particular  figure  of  merit  (such  as  the  A,  discussed  below).  In  either  of  the  above  conventions,  this 
reference  composition  becomes  a  function  of  film  thickness,  approaching  the  (stable  or 
metastable)  liquidus  composition  (Aj)  as  h  +oo.  An  alternative,  and  perhaps  theoretically  more 
elegant,  strategy  is  to  select  a  constant  reference  film  composition  (independent  of  h).  In  this  case, 
this  reference  composition  must  be  XL  to  ensure  that  crntctfacial  (+oo)  =  0 ,  even  if  the  average  film 
composition  can  be  significantly  different  from  XL  when  the  film  is  thin. 

Furthermore,  the  equilibrium  yc\  should  consider  the  effects  of  the  near-interface  gradients  in 
composition  (adsorption)  and  structure  (partial  ordering),  which  are  generally  difficult  to  quantify. 
An  useful  strategy  is  to  define  a  new  “un-relaxed”  yj0)  (>  yc\)  for  a  hypothesized  “step”  interface 
between  a  crystal  and  a  perfect  liquid  (without  any  adsoiption  and  near-interface  ordering),  which 
can  be  easier  to  quantify  (as  discussed  in  the  following  section).  Then,  we  can  define  a  new 
thermodynamic  variable  with  respect  to  a  hypothesized  uniform  film  of  a  perfect  liquid,  as: 
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(B.10) 


=  Y^~2Yf{Xfllm) 

amor film  ) 


Two  particular  Vs  that  can  be  quantified  are: 

<  (®<Xjilm  <1) 

A5W 


It  is  easy  to  demonstrate  that: 
A^  <  At  A,  <  A, , 

and 


lim  —  — >  1 

/l— >+oo  ^ 


a;  _  ^-2yri  • 
.4 


(B.l  1) 
(B.12) 


(B.l  3) 


All  three  A’s  represent  the  thermodynamic  tendency  for  stabilizing  subsolidus  liquid-like  IGFs, 
and  they  scale  the  actual  film  thickness.  Since  the  exact  form  of  the  interfacial  potential  is 
generally  unknown,  none  of  them  can  guarantee  more  prediction  power  than  others.  While  A* 
appears  to  be  conceptually  more  rigorous,  A  and  /w  are  generally  easier  to  quantify.  In  an  earlier 
publication  [6],  we  quantified  A  via  a  numerical  method.  More  recently,  we  derived  an  analytical 
solution  for  A,  for  regular  solutions  [7]  (and  in  this  appendix  we  further  generalize  it  to 
subregular  solutions),  which  is  practically  useful. 

In  summary,  A,  ,  which  is  the  maximum  thickness  of  a  stable  IGF  assuming  an  uniform  film 
composition  of  XL  and  no  interfacial  interactions,  appears  to  be  the  most  convenient 
thermodynamic  variable  to  be  used  practically.  Thus,  it  will  be  used  here.  To  compute  XL  as  a 
function  of  temperature  and  bulk  composition,  interfacial  energies  (y ’s)  and  bulk  free  energies  (to 
compute  AGamorph)  need  to  be  quantified  via  statistical  models  or  computational  thermodynamic 
methods,  which  are  discussed  in  following  sections. 


(B.3)  Estimation  of  Interfacial  Energies  via  Statistical  Models 

The  interfacial  energies  (y’s)  can  be  estimated  by  either  lattice-gas  models  or  Miedema  type 
“macroscopic  atom”  models. 


Lattice-Gas  Models: 


In  the  lattice  gas-model,  the  crystal-liquid  interfacial  energy  can  be  expressed  as  [8]: 


Td°')  -  Y<a>-{A]Xa  +  Y<b>-{b}xb  + 


vi1oL(XLAXSB+XSAXLB) 

V2/-i 


m  i<i>iXAXB  mi<i)sXAXB 


l/2/3 


1/2/3 


(B.l  4) 
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where  y<4>_{A)  and  y<B>_{B)  represent  the  solid-liquid  interfacial  energies  of  pure  A  or  B 

respectively;  XA  and  X'A  are  atomic  fractions  of  A  in  solid  and  liquid,  respectively;  similarly, 

XSB  and  Xlb  are  atomic  fractions  of  B  in  solid  and  liquid,  respectively  (  XSA  +  XB  =  1 ;  XLA  +  XB  = 
1);  coL  and  cos  are  the  regular  solution  parameters  for  solid  and  liquid,  respectively;  m{  is  the 
fraction  of  the  (liquid  type)  bonds  that  are  “cross”  the  interface;  V  is  molar  volume;  V213  is  molar 
area;  R  is  gas  constant;  and  T  is  temperature.  Essentially,  the  estimation  of  solid-liquid  interfacial 
energy  by  Eq.  (B.14)  is  through  calculating  the  bond  energies.  Elere  the  key  assumptions 
(simplifications)  are: 

•  The  solid-liquid  interface  can  be  represented  by  a  step  function  (i.e.,  adsorption  at  the 
interface  and  near-interface  ordering  are  not  considered);  thus  only  the  bonds  between  the 
first  layer  of  the  liquid  and  first  layer  of  the  solid  contribute  to  the  interfacial  energy; 

•  These  bonds  at  the  interface  is  liquid  type;  and 

•  The  reference  state  is  bulk  binary  solutions  (assuming  that  the  bulk  solid  is  in  a 
thermodynamic  equilibrium  with  the  bulk  liquid  phase). 


-B-B-  A-B-B-B-A-  B-  B-B-B 

- + — I — 4 — I — I — f — | — -I- — | — | — 4-- 

=A=  A=  A=A=A=A=A=A=  A=A=  A 

II  II  II  II  II  II  II  II  II  II  II 


>-  Liquid 
Interface 
>-  Solid 


Figure  B.l  Schematic  illustration  of  the  lattice-gas  model  for  a  liquid-solid  interface  of  an  A-B 
binary  system.  Assuming  that  the  solubility  of  B  in  solid  is  negligible  (XB  =  0). 


If  the  solubility  of  B  in  A  crystal  is  negligible  (XB  =  0;  XSA=I),  Eq.  (B.14)  can  be  simplified 
to: 


Yd  ~  Y<A>-{A) 


,  rnx(oL 


V 


2/3 


Kf 


(B.l  5) 


where  the  first  term  is  due  to  the  enthalpic  contribution;  the  second  term  is  the  interaction 
contribution  (a  “chemistry”  term  resulted  from  the  formation  of  a  different  number  of  A-B  bonds 
at  the  interfaces  as  compared  that  in  the  bulk  liquid).  The  above  lattice-gas  model  is  schematically 
illustrated  in  Fig.  B.  1 . 

The  above  model  can  be  further  modified  by  considering  the  adsorption  and  compositional 
gradients  at  the  solid-liquid  interfaces,  which  are  discussed  by  Shimizu  and  Takei  [9,10].  An 
entropic  term  associated  with  ordering  at  the  interfaces  (discussed  below  in  the  Benedictus  model) 
can  also  be  added. 
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“Macroscopic  Atom”  Models: 


Instead  of  using  the  lattice-gas  type  models  described  above,  we  adopt  and  modify  a 
Miedema  type  “macroscopic  atom”  model  to  estimate  interfacial  energies.  The  “macroscopic 
atom”  model  is  developed  on  a  similar  base  as  lattice-gas  models,  and  it  is  more  realistic  to 
represent  actual  transition  metal  alloys.  This  “macroscopic  atom”  model  has  the  following 
features  (that  are  useful  for  achieving  our  objectives): 

•  It  considers  the  different  molar  volumes  for  A  and  B; 

•  It  represents  average  interfaces  without  anisotropic  effects  (while  the  lattice-gas  models 
considers  the  anisotropic  effects  via  the  coordination  numbers);  and 

•  It  provides  a  systematical  method  to  evaluate  all  required  thermodynamic  parameters  (for 
all  binary  transition  metal  alloys  as  well  as  some  other  alloys). 

Benedictus  at  al.  reported  a  systematical  approach  to  evaluate  interfacial  energies  based  on 
the  Miedema  type  “macroscopic  atom”  model,  and  the  entropic  contribution  is  added  into  the 
model  [1 1],  In  this  model,  the  interfacial  energy  for  an  A  crystal  (assuming  that  the  solubility  of 
B  in  A  crystal  is  negligible)  and  a  binary  A-B  liquid  is  expressed  as: 


rrfuse 
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(B.16a) 


where  //(""’is  the  fusion  enthalpy  of  A;  Aft  jl,'/l«cc  's  the  enthalpy  of  solution  (similar  as  shown 
in  the  Miedema  model);  C0  ~  4.5  xlO8;  V  is  the  molar  volume  (neglecting  thermal  expansion);  R 
is  the  gas  constant;  and  Fs4  represents  the  “area”  fraction  A-B  bonds  at  the  interface,  which  is 
expressed  as: 


F  = 

1  B 


xbvb 


2/3 


(l-XB)VA2/3+XBV2/3 


(B.17) 


where  XB  is  the  atomic  fraction  of  B  in  the  liquid  phase. 

In  Eq.  (B.16a),  liquid-crystal  interfacial  energy  is  comprised  of  enthalpic,  interaction,  and 
entropic  contributions,  which  are  articulated  in  Ref.  [11]. 

It  is  important  to  note  that  the  interaction  term  in  Eq.  (B.16a)  is  derived  assuming  that  the 
reference  states  are  solid  pure  A  and  pure  B  (which  is  valid  for  the  discussion  of  solid-state 
amorphization  —  the  initial  objective  for  Benedictus  et  al.' s  work).  In  our  case,  however,  we  wish 
to  use  the  thermodynamic  equilibrium  state  as  the  reference  state.  To  achieve  this,  we  may 
assume  that  the  reference  state  is  set  by  the  bulk  liquid  (by  assuming  the  bulk  solid  A,  with 
negligible  solutes  of  B,  is  in  a  thermodynamic  equilibrium  with  the  bulk  liquid).  In  other  words, 
the  reference  state  for  the  interaction  term  is  set  by  a  hypothetical  liquid-liquid  interface,  which 
should  have  zero  excess  interfacial  energy  but  the  “macroscopic  atom”  model  would  produce  the 
following  interfacial  energy: 
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where  [2 F£  •  ( 1  —  F'/  )]  represents  the  “area”  fraction  of  A-B  bonds  for  this  hypothetical  liquid- 
liquid  interface  (which  should  produce  zero  excess  interaction  energy  by  definition).  Thus,  Eq. 
(B.16a)  should  be  revised  to 


rj- fuse 
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(B.16c) 


Or: 
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(B.16d) 


Here,  we  use  Mo-Ni  as  an  example.  Specific  to  this  system  of  our  interest,  the  liquid-crystal 
interfacial  energy  was  derived  from  Eq.  (B.16d)  and  determined  as 


ym-M°  =[0  138_0 .160(77^)2  +8.92x10  5r]J/m2 
where 
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(B.  19) 


The  average  GB  energy  of  Mo,  y^b ,  in  polycrystalline  materials  was  determined  by  using 
the  Turnbull’s  estimation  [12,13,14,15]: 

r™*i/3 -rL  (b.20) 

The  surface  energy  of  Mo  at  T=  OK  was  estimated  by  enthalpy  of  evaporation: 
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In  reality,  surface  energy  is  temperature-dependent,  where  both  the  thermal  expansion  effect  and 
entropy  effect  need  to  be  considered. 


rL(Tl)  = 


(jXTZ  +bMJl 

CO™1 


(B.22) 


bMo  is  a  materials-dependent  constant.  Furthermore,  structure  relaxation  at  free  surface  will  also 
reduce  its  excess  energy.  Nevertheless,  these  effects  are  relatively  small  (as  compared  with  the 
approximations  associated  with  Eq.  (B.20)),  and  they  are  generally  more  difficult  to  quantify. 
Thus  they  were  neglected  in  our  calculation  for  simplicity.  Combining  Eq.  (B.20)  and  (B.21) 

gives  an  estimation  of  y^b  «  1.0  J/nr. 


The  limitations  of  this  model  include: 

•  Adsorption  at  liquid-crystal  interfaces  is  not  considered.  This  is  addressed,  in  part,  in 
lattice-gas  model  by  Shimizu  &  Takei.[9,10] 
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•  Anisotropic  effects  are  not  considered.  This  is  considered  in  a  sophisticated  lattice-gas 
model  by  Wynblatt  et  al.  [16], 

•  Segregation  within  the  lattice  at  GBs  is  not  considered.  This  can  be  considered  by  Wynblatt 
et  al.’ s  model  [16], 

•  The  “asymmetric  effects,”  which  are  not  addressed  here  in  regular  solution  type  models,  is 
discussed  by  in  subregular  solution  model  by  Antion  and  Chatain  [17], 


(B.4)  Estimation  of  the  Free  Energy  Penalty  for  Forming  An  Undercooled  Liquid 

CalPhaD  Methods: 

The  Calculation  of  Phase  Diagram  (CalPhaD)  methods  were  adopted  to  determine  the  free 
energy  penalty  to  form  an  undercooled  liquid  (AGam0rph).  CalPhaD  is  a  well  established  method  to 
compute  bulk  phase  diagrams  from  (usually  empirically -fitted)  thermodynamic  functions  of  free 
energies.  The  Gibbs  free  energy  of  a  phase  0  in  a  binary  A-B  system  can  be  expressed  as: 

G®  =  ^X,.-°G®  +RT  ln^+^Ga>  (B.23) 

i=A,B  i=A,B 

where  °G,®  is  the  Gibbs  free  energy  of  the  pure  element  i  (=  A  or  B)  that  is  present  in  0  phase;  X, 

is  the  atomic  fraction  of  element  i;  xs  G'"  is  the  excess  Gibbs  free  energy,  which  can  be 
empirically  expressed  in  a  Redlich-Kister  polynomial: 

(B.24) 

7=0 

Here,  L®  ’s  are  empirically-fitted  parameters  for  the  phase  0.  If  n  =  0,  the  phase  0  is  a  regular 
solution,  and  Eq.  (B.24)  is  reduced  to 

■KG'"  =  L%XaXb  =  coXAXB  .  (B.25) 

where  co  (=  L® )  is  the  regular  solution  parameter.  If  co  =  0,  the  phase  0  is  an  ideal  solution 
(®G®  =  0). 

If  n  =  1,  the  phase  0  is  a  sub-regular  solution.  Eq.  (B.24)  can  be  rewritten  as 

=  L%XaXb  +LfXAXB(XA  -AY).  (B.26) 

Gibbs  free  energy  functions  for  compounds  and  ordered  solutions  can  be  constructed  using 
different  models,  which  are  discussed  in  Ref.  [18].  For  a  given  binary  A-B  system,  the  Gibbs  free 
energy  functions  can  be  developed  for  all  possible  phases.  Then  a  bulk  phase  diagram  can  be 
constructed  by  minimizing  the  total  free  energy  of  the  system.  In  two -phase  regions,  the 
equilibrium  states  can  be  found  graphically  via  a  well-established  “common  tangent  construction” 
method. 

For  the  Ni-Mo  binary  system  (as  the  first  computation  example),  four  phases  (BCC,  FCC,  8- 
NiMo,  and  liquid  phases)  are  considered.  The  excess  Gibbs  energies  of  them  were  determined 
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and  the  Gibbs  free  energies  of  different  phases  are  listed  in  Table  B.l(a)  [19]  and  Table  B.l(b) 

[20]. 


Table  B.l(a)  Gibbs  free  energies  of  the  different  phases  in  the  Ni-Mo  system  [19]. 


Phase 

Gibbs  free  energy,  J/mol 

Liquid  phase 

+Xmo°G'1  -TASconf  +  XNiXMo  [-46540  - 1 9.53T  -  29 1 5(XMo -XNi)] 

Mo-rich  BCC 

XNi  °GBfc  +  XMo  °Gbmc0c  -  TASconf  +  46422X NiX Mo 

Ni-rich  FCC 

XNi°G%c+XMo°G™c  - TASconf  +XNiXMo[-4m.7-5.96T-mmXMo -XNi)] 

6-NiMo 

24° GmC  +  32 ”G*C  -  2 1 2 1 00  + 1 089J  - 1 42 J ln(T) 

Note.  A>S'con p  XW]lnXK)  X\]()lnX\]<r 


Table  B.l(b)  Gibbs  free  energies  of  the  different  phases  in  the  Ni-Mo  system  [20]. 


Phase 

Gibbs  free  energy,  J/mol 

Liquid  phase 

XNi°G'“‘+XMo°G,*0-TASconf  + 

XNiXMo[-39597  +  \5.935T-7373(XMo-XNi)  +  (-m23  +  5.55lT)(XMo-XNi)2] 

Mo-rich  BCC 

XNi°GBNfc  +XMo°GZc -TAScon/:  +  XNjXMo [2769 1  +  1 8792(XMo  -XK)] 

Ni-rich  FCC 

y  o  s-iFCC  .  y  orFCC  rp  A  o  ,  y  y 

A  Ni  Ni  “l_AMo  ^Mo  1  ^conf.  +  A  Mo  * 

[-8916  +  3.591T  +  (5469  - 0.249T)(XUo  - XNi)  +  (-1549 - 2.14\T){XMo  - XNi)2] 

S-NiMo 

24°G^C  +  20°G^oc  +  12°  G^f  - 169981 +  1154.981T- 155. 4847’ln(T) 

+  24° GBfc  +20° G^c  +12° G^cf  -154106  +  2855.00 IT -394.923T\n(T) 

-199.856 T  +  0.3367  •  [(-829211  +  825.923 T)  +  0.3571  •  (-417368  +  326.504T)] 

Note.  AScon p  X^\\hiX^\X  Xy[0lnXy[0. 


Table  B.l(c)  Comparison  of  phase  diagram  data  comparison  by  the  two  groups  of 
thermodynamic  functions  listed  in  Table  B.l(a)  and  Table  B.l(b),  respectively. 


Key  markers 

Experimental  value 

Krisk’s  [19] 

Zhou  et  al.’s  [20] 

T  °C 

peritectic? 

1362  [21,22] 

1364 

1345 

Solidus  line  and  solvus 
line  of  BCC 

From  ref.  [23,24] 

Good  fit 

No  validation 

The  thermodynamic  functions  in  Table  B.l(b)  is  more  sophisticated  than  those  in  Table  B.l(a) 
by  considering  entropic  contribution  to  the  Gibbs  free  energy  of  8-NiMo  compound. 

Consequently,  the  thermodynamic  functions  listed  Table  B.l(b)  achieved  better  accuracy  for 
predicting  FCC  phase  and  intermetallic  compounds.  However,  the  thermodynamic  functions  in 
Table  B.l(a)  lead  to  better  predictions  for  the  solidus  and  solvus  lines  of  the  Mo-rich  BCC  and 
the  peritectic  temperature,  which  are  our  primary  concerns.  Therefore,  we  used  Table  B.  1  for  our 
calculations. 
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When  the  thermodynamic  functions  of  all  phases  are  known,  we  can  construct  binary  phase 
diagram  by  drawing  common  tangent  line  and  finding  equilibrium  phases.  Constructions  for  two 
representative  temperatures  of  the  Mo-Ni  binary  system  are  shown  in  Fig.  B.2. 


X(Ni)  X(Ni) 

(b) 


Figure  B.2  Schematic  illustration  of  phase  diagram  constructions:  (a)  The  Mo-Ni  binary  phase 
diagram,  (b)  Plots  of  Gibbs  free  energies  of  different  phases  vs.  Ni  atomic  fraction  for  1200  °C  (< 
Tperitectic)  and  1400  °C(>  Tperitectic),  in  which  common  tangent  lines  were  drawn  to  find  the  phase 

boundaries. 
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Determination  of  AGamorph  and  A,L: 

The  CalPhaD  method  was  used  to  compute  AGamorph  and  then  /..  Again,  we  used  Mo-Ni  as 
the  fist  computation  example.  From  now  on,  the  subscript  “Ni”  is  dropped  for  brevity,  all  the 
X’s  (X,  Xflim,  Xs,  Xj)  in  equations  are  referred  to  as  atomic  fraction  ofNi  in  Ni-Mo  binary  system 
with  the  following  definitions. 

X:  the  Ni  fraction  (general) 

Xfllm:  the  Ni  fraction  of  IGF 

X0 :  the  Ni  fraction  of  the  bulk  (BCC)  phase  (which  sets  the  bulk  chemical  potentials) 

Xs:  the  Ni  fraction  on  the  solidus  line  or  (if  T  <  Tv,cllca,c)  its  metastable  extension 

XL:  the  Ni  fraction  on  the  liquidus  line  or  (if  T  <  Tpm{cclK)  its  metastable  extension 

Xsolvus:  the  Ni  fraction  on  the  solvus  line  (for  T  < 

In  particular,  analytical  expressions  for  AL  can  be  obtained  for  Ni-doped  Mo  with  two 
assumptions.  First,  the  IGF  has  a  uniform  composition  that  is  equal  to  the  composition  on  the 
liquidus  line  or  its  metastable  extension.  Second,  the  solid  (BCC)  phase  is  a  regular  or  sub¬ 
regular  solution  so  its  mixing  Gibbs  energy  is  given  by 


A GZid  =  RT\X0  hi  X0+(\-X0)  ln(l  -  X0 )]  +  L%X0 (1  -X0)  +  L> X0 (1  -  X0 )(1  -2X0), 


(B.27) 


where  X0  is  the  Ni  atomic  fraction  in  the  solid  phase;  Ls0  (=co)  and  L\  are  sub-regular  (or  regular) 

solution  parameters.  In  our  calculations,  the  solid  phase  is  the  BCC  phase  because  the  FCC 
phase  is  unstable  in  this  composition  region  of  interest. 

First,  we  will  determine  the  AGamorph  in  the  single-phase  (BCC)  regime.  Fig.  B.3  shows  an 
example  at  T=  1400  °C  (T>  TpaiU:clK),  in  which  we  assume  that  the  composition  of  the  bulk  BCC 
phase  (X0)  is  lower  than  the  bulk  solidus  composition  (Xs).  The  AGamorph  is  defined  as  the 
difference  between  formation  energy  of  the  liquid  phase  (the  red  line  in  Fig.  B.3(a))  and  the 
reference  state  set  by  the  chemical  potential  of  the  bulk  phase  of  a  composition  X0  (the  black 
straight  line  labeled  in  Fig.  B.3(a)).  This  reference  state  is  given  by: 


GRef(X)  =  X 


film 


■Mm  +(WXfJ.Muo  =GfbJX„)  +  (Xfllm -X„) 


dGL 


dX 


x=x„ 


which  corresponds  to  the  tangent  line  extended  from  X0  in  Fig.  B.3  (a).  Thus,  the  free  energy 
penalty  to  form  a  metastable  liquid  of  composition  XL  (being  labeled  by  a  purple  double  arrow  in 
Fig.  B.3(a))  is  given  by: 


'liq' 


GL(X o  +  (xl-x„) 


dGL 


dX 


X=xa 


(B.29) 


where  Gjiq 


and  G/cc  are  the  formation  free  energies  of  liquid  and  solid  BCC  phases,  respectively. 
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Gibbs  free  energy,  J/mol 


(a) 


(b) 

Figure  B.3  (a)  The  representation  of  AG'amoiph  (XL)  in  the  plots  of  Gibbs  free  energies  vs.  Ni 

fraction  for  a  case  of  X0  <XS  and  T=  1400  °C  (  >  Tperitectic  ),  where  X0  is  the  bulk  composition.  Xs 
and  XL  are  equilibrium  solidus  and  liquidus  compositions,  respectively,  (b)  The  corresponding 
bulk  phase  diagram,  where  the  single-phase  regime  of  interest  is  shaded. 


By  definition,  AGamorph  vanishes  if  the  composition  of  solid  BCC  phase  is  on  the  bulk 


53 


solidus  line  (X0  =  XS),  which  is  the  case  represented  the  common  tangent  line  between  Xs  and  XL 
in  Fig.  B.3  (a);  thus 


(%(*,)  =  G'k(Xs)  +  (Xl-Xs) 


dG, 


bcc 


dX 


(B.30) 


Combining  Eqs.  (B.27-B.30)  while  assuming Xmm  =  X,  gives: 


AG amorph  =  RT 


X.in^  +  a-XJin^ 

Aq  I  A  Q 


As  - X0)[2 Xl(Lb0cc  +  3 LBCC) - (Lb0cc  +  3 L°cc  +  6 L°CC){XS  +  X0)  +  4Lfcc(X2  +  XSX0  +  X20)\ 

(B.3 1 ) 


If  the  solid  BCC  phase  is  a  regular  solution  ( LBCC  =  0;  LBCC  =co),  it  can  be  reduced  to 


AG  ,=RT 

amorph 


XL  In  +  (1  -  XL  )  In  1  X 


Xn 


1-Xn 


co(Xs-X0)(  2XL-XS 


■x0). 

(B.32) 


Second,  in  the  single  BCC  phase  regime  below  rperitectic,  the  calculation  is  almost  the  same  as 
that  described  above  for  T  >  Tpa itectic-  The  same  equations  are  used,  and  the  only  difference  is  that 
Xs  and  XL  are  now  the  Ni  fractions  on  the  metasatble  extensions  of  the  solidus  and  liquidus  lines, 
as  illustrated  in  Fig.  B.4  (for  an  example  for  T=  1200  °C).  Note  that  now  the  solid  solubility 
corresponds  to  the  solvus  line  2fso|vus  (instead  of  As,  which  is  on  the  metastable  extension  of 
solidus  line)  for  T>  Tperitectic. 


Finally,  in  the  sub-peritectic  two-phase  regime,  AGamorph  (XL)  is  a  constant  at  a  specific 

temperature  and  it  does  not  depend  on  the  overall  bulk  composition  (because  the  bulk  chemical 
potential  is  a  constant.  Therefore,  AGamorph  (XL)  can  be  determined  by  using  Eq.  (B.32)  or  (B.31) 


and  assuming  X0  =  Aso|vus: 


AG  amorph  =  RT 


XL  In  f-  (1  -  X  L  )  In :  1  ^ 


X., 


\-X, 


for  a  regular  solution,  or 


~oAXs  -  X solvus)(2X L  -Xs  -  XsolvJ 

(B.33) 


AG amorph  =  RT 


XL  In  r  (1  - XL)  In  ~~~~ 

^  solvus  v  solvus 


As  -  Xsolmi2 XL{Lfc  +  3 Lfcc)  -  (. Llcc  +  3 Lbcc  +  6 Lfcc)(Xs  +  Xsolvus)  +  4Lfcc(X2  +  XSX  +  X2ohj\ 

(B.3  4) 

for  a  sub-regular  solution. 


Furthermore,  in  the  sub-peritectic  two-phase  regime,  TL’s  are  horizontal  lines,  and  lines  of 
constant  AL  can  be  constructed  graphically  by  finding  the  intersections  of  constant  AL  lines  and 
the  solvus  line,  as  illustrated  in  Fig.  B.5(c). 
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-50000 


Afsolvus 


(a) 


(b)  (c) 

Figure  B.4  (a)  The  representation  of  AGamorph  (X\_)  in  the  plots  of  Gibbs  free  energies  vs.  Ni 

atomic  fraction  for  a  case  of  X0  <  Aso|vus  and  T=  1200  °C  (  <  rperi,ectic),  where  Xs  is  the  composition 
on  the  metastable  extension  of  the  solidus  line;  (b)  and  (c)  are  the  corresponding  bulk  phase 
diagram  and  an  expanded  view.  Bulk  solidus  and  liquidus  composition  (Xs  and  XL  respectively) 
are  extended  into  metastable  region  as  dashed  line. 
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Gibbs  free  energy,  J/mol 


-50000 


(a) 


(b)  (c) 

Figure  B.5  (a)  Illustration  of  determining  AGamorph  at  T=  1200  °C  in  the  sub-peritectic  two 

phase  regime;  (b)  Ni-Mo  binary  phase  diagram.  The  red  dashed  line  represents  regime  of 
calculation,  (c)  Construction  of  lines  of  constant  /lL. 
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(B.5)  Computation  Resulted  for  Ni-Doped  Mo 


We  plotted  lines  of  constant  computed  AL  in  the  bulk  phase  diagram  of  the  Mo-Ni  system. 

This  represents  a  computed  “GB  diagram,”  where  AL  represents  the  maximum  thickness  of  a 
quasi-liquid  IGF  of  composition  XL  that  can  be  stabilized  at  Mo  GB  (an  average  random  GB) 
without  the  consideration  of  interfacial  forces  and  through-thickness  compositional  and  structural 
gradients.  It  represents  the  thermodynamic  tendency  for  a  Mo  GB  to  disorder.  AL  (or  A)  is  not  the 
actual  IGF  thickness,  but  it  should  scale  the  thickness  of  the  actual  (thermodynamically  stable) 
quasi-liquid  IGF. 

Fig.  B.7  further  shows  computed  of  AL  and  A  vs.  Ni  atomic  fraction  at  three  representative 
temperatures.  The  same  trends  are  found  for  AL  and  A  (although  AL  <  1);  and  both  are  divergent  as 
bulk  Ni  composition  approach  to  the  solidus  line  (or  its  metastable  extension).  For  T<  Tpcnlcaic, 

AL  and  A  level  off  at  the  bulk  solvus  line  X0  =Xsoivus,  above  which  the  bulk  chemical  potential  is  a 
constant  (despite  any  change  in  the  overall  bulk  composition  in  the  BCC-NiMo  two  phase  region). 

In  the  single-phase  region,  computed  Ah  (or  A)  increases  with  increase  in  bulk  composition. 
Above  the  peritectic  temperature,  computed  AL  (or  A)  is  divergent  at  the  bulk  solidus  temperature. 
Below  the  peritectic  temperature,  computed  AL  (or  A)  levels  off  at  the  bulk  solvus  temperature,  and 
it  becomes  a  constant  in  the  BCC-NiMo  twp-phase  regime. 


Ni  atomic  fraction 

Figure  B.6  Computed  “grain  boundary  diagram”  for  Ni-doped  Mo,  where  computed  lines  of 
constant  AL  are  plotted  in  the  binary  bulk  phase  diagram.  Note  that  Ah  is  the  maximum  thickness 
of  a  quasi-liquid  IGF  of  composition  XL  that  can  be  stabilized  at  an  average  (random)  Mo  GB 
without  the  consideration  of  interfacial  forces,  and  it  represents  a  thermodynamic  tendency  for  an 

average,  random,  Mo  GB  to  disorder. 
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A 

AL 


(a) 


Ni  atomic  fraction 


(b) 


Ni  atomic  fraction  Ni  atomic  fraction 


(C) 

Figure  B.7  Plots  of  XL  and  X  vs.  Ni  atomic  fraction  at  (a)  1300  °C  (a),  (b)  1400  °C  and  (c)  1495  °C, 
respectively.  The  dashed  vertical  line  in  (a)  represents  bulk  solvus  composition;  the  solid  vertical 
lines  represent  the  bulk  solidus  compositions  (or  the  composition  on  the  metastable  extension  of 

the  solidus  line  for  1300  °C  in  (a)). 
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(B.5)  Computation  Results  for  Several  Doped  W  Alloys 


Similar  computations  can  be  conducted  for  other  binary  transition  metal  systems.  Computed  “GB 
diagrams”  (noting  that  these  are  “A  diagrams”)  of  Pd,  Ni,  Fe,  and  Co  doped  W  are  shown  in  Fig.  B-8. 
These  computations  were  conducted  using  the  CalPhaD  data  listed  in  Table  B.2.  Fig.  2  and  Fig.  9  in  the 
narrative  are  colored  versions  of  this  GB  diagrams,  respectively.  The  correctness  of  these  computations  is 
validated  by  comparing  with  experimentally  observed  onset  sintering  temperatures  and  direct  FIRTEM 
characterization,  which  are  discussed  in  the  narrative  of  this  technical  report.  Computation  has  also  been 
conducted  for  Cu-doped  W;  the  results  shown  that  IGFs  cannot  form.  Consistently,  it  was  well  known  that 
Cu  has  no  effect  for  activated  sintering  of  W. 
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Figure  B.8:  Computed  “GB  diagrams”  (A  diagrams)  of  (a)  Pd-doped  W,  (b)  Ni-doped  W,  (c)  Fe-doped  W, 
and  (d)  Co-doped  W.  In  each  bulk  phase  diagram,  a  series  of  dashed  and  dotted  lines  were  plotted  to 
represent  the  conditions  that  4,  2,  1  and  0.5  nm  thick  quasi-liquid  IGFs  can  be  thermodynamically  stable  at 
GBs.  The  ranges  of  subsolidus  activated  sintering  are  also  labeled. 

Reprinted  from  an  AFOSR  supported  publication  [Curr.  Opin.  Solid  State  Mater.  Sci.  12:81-88  (2008)] 
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Table  B.2  Binary  thermodynamic  functions  used  in  the  calculations. 


Phase 

Gibbs  free  energy  (J-mol1) 

Ref. 

Liquid 

XNLG^+XW°G^  -  TASconf  +  XNjXw[l  6290  - 1 0.257  -  22450(XM -Xw)] 

BCC 

XNi°GbN?  +  Xw°Gb™  - TASconf  +XNiXw(lM12) 

W-Ni 

FCC 

Xm°G%  +XW°G^  —  TASconf  +XNiXw[2556  +  U.6T-52900(XNt-Xw)]  +  Gmag. 

25 

jNi4W 

fGf™  +\°Gb‘c  -  2572.4  -  0.94947 

j  NiW 

\°Gf£  +YG^CC  -  4442.95  + 1 .57 

j  NiW2 

f  G%  +| °Gbvcc  - 2963.02  +  7 

Liquid 

XCo°GY  +  Xw°G%  -TASconf  +  XCoXw (-42606.6  + 1 5.50647') 

W-Co 

BCC 

YCo  °Gbcc;  +  Xw  "G'Y  -  TASconf  +  XCoXw  (54648) 

26 

FCC 

XCo°GY  +  Xw°G*°  -TASconf  +  XCoXlv[-] 3038  +  8.0387 -  1 1 832(YCo -Xw)\ 

& Co6W7 

V3°Gc7  +i \°Ghwcc  -  3796.9  -  0.65267 

Liquid 

Xr°G%  +Xw°G'p  —  TASconj-  +  XFeXw\—2>691  +  5.657  —  5450(XFc -Xw)] 

BCC 

XFe°Gt  +  Xw"GhY  - TASconf  +  X FeXw\4l544 -  12621(XFe  - Xw)]  +  Gmag_ 

FCC 

XFe°Gf;:  +Xw°Gp  —  TASconf  +  XFeXw[25977  -  3.36T-7258.6(XFe - Xw)]  +  G,mg. 

W-Fe 

|i-phase 

y Fe  GFe..W:Fe  +  JV  GFeWW  +  ^(Ppe  ^  47%  +  JV  ^  JV) 

where:  yF  and  y^,  are  the  site  fraction  of  element  Fe  and  W  on  sublattice  3. 

°GMFeW:Fe  =  V°Gg  +2°GYc  +4°Gy“7  -14300  +  26.77)/13 

27, 

28 

°Gl„.w  =  (TGg  +2 °Gbwcc  +4 °GbY  -  53450  + 1 97)  / 1 3 

j  Fe2W 

\°G%  +\°GbY  -  4500  + 1 .666677 

Liquid 

XPd°G%  +  XW°G»  - TAS^  +  XPdXw (13868.7  + 10.93227) 

W-Pd 

BCC 

Xm°G%+Xw°G%°  -TASnmf  +  XPrfX(r(9386.8  + 14.63297) 

29 

FCC 

XPd°G%  +XW°GY -TASconf  +  XPdX„(-13144  +  16.68667) 
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Appendix  C: 

Sintering  and  GB  Diffusivities  of  Ni-Doped  Mo 


This  appendix  documents  a  study  to  investigate  the  sintering  behaviors  of  Ni-doped  Mo,  in 
which  we  extracted  GB  diffusivities  systematically.  The  objectives  are  twofold: 

•  We  extracted  GB  diffusivities  to  provide  indirect,  yet  systematical,  experimental 
evidence  to  systematically  validate  the  thermodynamic  models  and  computed  GB 
diagrams  developed  in  this  AFOSR-Y1P. 

•  Good  model-experimental  agreements  provided  further  evidence  to  quantitatively  support 
the  proposed  solid-state  activated  sintering  mechanism,  i.e.,  the  subsolidus  accelerated 
sintering  is  due  to  the  enhanced  diffusion  in  quasi-liquid  IGFs. 

The  experiment  al  procedures  are  described  in  Appendix  A.  The  XRD  confirmed  that  all  the 
Ni  was  transformed  to  S-NiMo  compound  (Fig.  C.l),  excluding  the  possibility  of  forming  a 
transient  liquid. 

Sintering  temperatures  were  1220  °C,  1300  °C,  1350  °C,  1400  °C,  1450  °C,  and  1495  °C, 
respectively,  while  the  peritectic  temperature  is  1362  °C.  Different  Ni  atomic  fractions  were 
selected  so  that  it  covered  all  three  phase  regimes.  The  selected  experimental  conditions  are 
schematically  illustrated  in  Fig.  C.2. 
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Figure  C.l  XRD  spectra  of  Ni-doped  Mo.  Almost  all  the  Ni  has  been  transformed  to  S-NiMo 
compound  after  presintering  at  1 000  °C  for  1  hour. 
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Ni  atomic  fraction 


Figure  C.2  The  selected  sintering  temperatures  and  compositions  are  marked  in  the  Mo-Ni  binary 

phase  diagram. 


(C.l)  Densification  and  Grain  Growth 


The  densification  data  for  isothermal  sintering  at  eight  temperatures  are  documented  in  “§C.5 
Tables  of  Densification  Data”.  Several  representative  plots  of  relative  density  vs.  sintering  time 
are  shown  in  Fig.  C.3  and  Fig.  C.4.  The  densification  rate  generally  increases  with  increasing  Ni 
content  or  temperature. 

The  error  bars  were  determined  by  following  methods.  We  selected  a  few  representative 
conditions  to  prepare  3-4  specimens  to  obtain  means  and  standard  deviations  of  the  measured 
densities  (to  represent  specimen-to-specimen  variations).  Based  on  these  measurements,  we 
found  that  the  standard  deviation  is  about  0.005  of  the  relative  density  (or  0.0013  in  the  linear 
shrinkage,  A L/L0)  for  specimens  generally.  These  error  bars  were  then  used  as  estimated  errors 
bars  for  other  measurements. 

The  measured  grain  sizes  of  samples  with  different  Ni  contents  and  sintering  temperatures  are 
shown  in  Table  C.2.  The  grain  growth  is  almost  negligible  (within  the  range  of  experimental 
errors)  for  moderate  doping  level  (<  1.5  at  %  Ni)  or  at  low  temperatures  (<  1300  °C).  For  the 
regime  where  the  grain  growth  is  not  negligible  (the  region  that  is  the  left-bottom  comer  by  the 
red  line  in  Table  C.2),  the  effects  of  grain  growth  were  considered  when  we  extracted  the  GB 
diffusivity  data  via  fitting  the  sintering  models. 
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1495  °C 


4%Ni 

- 2%Ni 

h— 1.7596NI 
-#-1.5%Ni 
l%Ni 
0.75%Ni 
— *-0.5%Ni 
—♦—Pure  Mo 
— *-0.25%Ni 


— *-4%Ni 
2%Ni 

- 1.75%Ni 

- 1.596NI 

— l%Ni 
-H—  0.75%Ni 
— *-0.5%Ni 
— 0.375%Ni 
-■-0.25%Ni 
♦  Pure  Mo 


(a) 


(b) 


Figure  C.3  Relative  density  vs.  sintering  time  at  (a)  1490  °C  and  (b)  1400  °C  (rpreitectic  =  1362  °C). 


4%Ni 

- 2%Ni 

1.25%Ni 
-#-l%Ni 
0.75%Ni 
-*-0.5%Ni 
— 0.375%Ni 
-«-0.25%Ni 
♦  Pure  Mo 


1300  °C 


1.7596NI 

0.75%Ni 

4%Ni 

0.5%Ni 

0.25%Ni 

Pure  Mo 


(a) 


(b) 


Figure  C.4  Relative  density  vs.  sintering  time  at  (a)  1350  °C  and  (b)  1300  °C  ( 7'prcitcctic  =  1362  °C). 


64 


Table  C.l  Measured  grain  sizes  of  Ni-doped  Mo  samples. 


(C.2)  Determination  of  the  Controlling  Densification  Mechanism 


The  GB  diffusivity  can  be  extracted  from  sintering  data  only  if  the  densification  is  controlled 
by  GB  diffusion.  We  can  determine  whether  densification  is  controlled  by  GB  diffusion  (vs.  bulk 
diffusion)  by  obtain  the  exponent  m  for  the  power  law  fitting  of  the  initial  stage  densification 
kinetics.  According  to  Johnson’s  model  [2]  and  Kang’s  sintering  theory  [3],  the  linear  shrinkage 
(A L/Lq)  in  the  initial  stage  follows  a  power  law: 


^L=  K  y  D-D 
4  [_  kTQP 


(C.l  a) 


The  experimentally  fitted  exponent  number  m  can  be  used  to  discriminate  what  the  sintering 
mechanism  is  in  control  [2].  If  m  =  -0.4-0. 5,  the  densification  is  controlled  by  bulk  diffusion, 
where  K  =  6,p  =  3>,  and  D  =  Dxl ■  Eq  (C.l a)  becomes 


^  Y ji/2 

4  kTG3  \ 


(C.lb) 
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If  m  =  -0.33,  the  densification  is  controlled  by  GB  diffusion,  where  K=  12,  p  =  4,  and  D  =  8 DGB. 
Eq.  (C.la)  becomes 

L0  [_  kTG4 

The  index  m  can  be  determined  experimentally  by  finding  the  slope  the  double  logarithmical  plot 
of  linear  shrinkage  vs.  time  via  linear  regressions,  which  are  illustrated  in  Fig.  C.5. 

The  exponents  (m’s)  were  calculated  and  listed  in  Table  C.2.  Only  in  two  samples  (pure  Mo 
and  0.25  at  %  Ni  doped  Mo,  both  sintered  at  1495  °C),  the  densification  was  controlled  by  bulk 
diffusion,  where  GB  diffusivity  were  not  measurable.  These  two  specimens  are  labeled  in  Fig. 
C.6  by  “X”.  In  all  other  specimens,  the  densification  was  GB  diffusion-controlled,  and  the  GB 
diffusivities  can  be  extracted.  This  sintering  mechanism  map  (Fig.  C.6)  is  consistent  with  a  prior 
report  [4]. 


.1/3 


(C.lc) 


X  1.5%Ni 
■  1%  N  i 

•  0.5%Ni 

♦  Pure  Mo 


Figure  C.5  Double  logarithmical  plots  of  linear  shrinkage  vs.  sintering  time  for  samples  sintered 
at  1495  °C.  The  slopes  correspond  to  the  exponents  (in's). 


Table  C.2  Fitted  exponents  (m’s) 


^^Ni% 

Index 

Pure  Mo 

0.25% 

0.5% 

1.0% 

1.5% 

1495  °C 

0.44  ±  0.06 

0.37  ±0.02 

0.23  ±  0.03 

0.32  ±0.02 

0.22  ±0.03 

1450  °C 

0.31  ±0.11 

0.23  ±  0.03 

0.22  ±0.03 

0.29  ±0.02 

0.20  ±0.03 

1400  °C 

0.29  ±  0.06 

0.24  ±0.01 

0.26  ±0.02 

0.25  ±  0.04 

0.23  ±  0.07 

1350 °C 

0.28  ±0.03 

0.23  ±  0.03 

0.31  ±0.02 

0.36  ±0.03 

1300 °C 

0.29  ±  0.06 

0.23  ±  0.02 

0.28  ±0.07 

1220  °C 

0.29  ±0.06 

0.35  ±0.03 

0.26  ±0.07 

0.33  ±0.08 
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Ni  atomic  fraction 


Figure  C.6  The  experimental  conditions  are  marked  in  the  Mo-Ni  bulk  phase  diagram.  The  blue 
“  indicates  where  densification  is  controlled  by  bulk  diffusion.  The  region  where  GB 
diffusivities  were  determined  by  fitting  the  initial  stage  sintering  model  and  the  region  where 
grain  growth  is  not  negligible  are  also  indicated  in  this  diagram. 


(C.3)  Extraction  of  GB  Diffusivities 


The  sintering  models  had  been  discussed  systematically  in  Appendix  A.  We  extracted  GB 
diffusivities  from  densification  data.  When  it  is  possible,  the  intermediate-stage  sintering  data 
(linear  shrinkage  >  4%)  are  used  for  better  accuracy.  In  the  case  of  pure  Mo  and  low  doping 
levels  (Ni  atomic  %  <  0.5),  the  total  densification  is  not  sufficient  to  apply  the  intermediate -stage 
sintering  model,  the  initial  stage  densification  data  were  used  to  obtain  GB  diffusivities. 

For  a  few  cases  where  both  initial  and  intermediate  stage  sintering  data  exist,  both  Johnson 
(initial  stage)  and  Coble  (intermediate  stage)  models  are  used  to  extract  GB  diffusivities  and  the 
results  are  listed  in  Table  C.3.  The  GB  diffusivities  determined  by  both  models  are  rather 
consistent. 

We  also  compared  the  GB  diffusivities  determined  by  the  German  model  (for  liquid-phase 
sintering)  and  the  Coble  model  (for  nominal  GB  diffusivity  via  solid-state  sintering  model)  in 
Table  C.4.  It  can  be  shown  that  these  two  models  were  roughly  equivalent.  The  calculated  results 
from  these  two  models  are  in  fact  consistent. 
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Table  C.3  Comparison  of  GB  diffusivities  extracted  from  the  Johnson  initial-stage  model  and  the 

Coble  intermediate -stage  sintering  model. 


T,°C 

1220 

1300 

1350 

1400 

1450 

1495 

Ni  at.  % 

0.25 

0.5 

0.5 

0.5 

0.75 

0.5 

From  initial 
stage  data 

1.89E-18 

0.97E-17 

4.81E-18 

7.13E-18 

1.74E-17 

9.27E-18 

From 

intermediate 
stage  data 

1.34E-18 

1.20E-17 

4.96E-18 

6.01E-18 

2.36E-17 

3.18E-18 

Table  C.4  Comparison  of  GB  diffusivities  determined  by  using  the  German  and  the  Coble  model 
for  specimens  sintered  in  the  liquid-phase  sintering  region. 


T,  °C 

1450 

1450 

1495 

1495 

Ni  at.  % 

2 

4 

2 

4 

The  German  model 

5.21E-16 

9.69E-16 

6.70E-16 

1.58E-15 

The  Coble  model 

9.84E-16 

13.7E-16 

9.86E-16 

1.21E-15 

(C.4)  Results  and  Comparisons 

Plots  of  extracted  GB  diffusivities  vs.  Ni  contents  for  six  different  temperatures  between 
1220  °C  and  1495  °C  are  shown  in  Fig.  C.7. 

Plots  of  extracted  GB  diffusivities  vs.  temperatures  are  shown  in  Fig.  C.8  (for  2  at.  %  Ni,  in 
the  two-phase  regions)  and  Fig.  C.9  (0-1  at.  %  Ni,  in  the  single-phase  region). 
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1495  °C 


Ni  atomic  % 


1450  °C: 


Ni  at.  % 

0.14 
0.12  - 
0.10  - 
3  0.08 
<  0.06 
0.04  - 
0.02 
0.00 
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12  3  4 

Ni  atomic  % 


1400  °C: 


Ni  atomic  % 


Ni  at.  % 


1350 °C: 
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1300 °C: 


1220  °C 


Figure  C.7  (Left)  Plots  of  GB  diffusivities  vs.  Ni  contents.  (Right)  Plots  of  linear  shrinkages 
(after  2-hour  isothermal  sintering)  vs  Ni  contents.  The  solid  vertical  lines  represent  bulk  solidus 
composition  (for  T>  Tpi itectic  =  1362  °C)  or  solvus  composition  (for  T <  7pritectic  =  1362  °C).  the 
dashed  vertical  lines  represent  metastable  solidus  composition  (for  T <  T]irmcUc  =  1362  °C). 


Temperature,  °C 


(a) 


Figure  C.8  (a)  GB  diffusivity  vs  temperature  for  2  at.  %  Ni  doped  Mo  (in  two-phase  regions: 
BCC  +  S-Ni-Mo  for  T <  rpritectic  or  BCC  +  liquid  for  T  >  Tpritectic). 

(b)  Associated  computed  1L  vs.  temperature. 
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0.5%Ni 

♦ 
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1400 
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0.03 
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Figure  C.9  Temperature  dependence  of  linear  shrinkage  (after  2  hrs)  and  GB  diffusivity  in  the 
single-phase  regime  for  (a)  pure  Mo,  (b)  Mo  +  0.25%  Ni,  (c)  Mo  +  0.5%  Ni,  and  (d)  Mo  +  1%  Ni. 
The  plots  of  2l’s  vs.  temperatures  are  also  shown. 
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The  comparisons  between  measured  GB  diffusivities  and  the  computed  GB  diagram  are 
summarized  in  Fig.  C.  10.  The  major  jumps  in  GB  diffusivities,  which  are  labeled  by  the  red  stars 
in  Fig.  C.  10(A)  and  (B),  occur  at  a  distinct  “GB  solidus  line”  that  corresponds  to  a  computed  A,L 
value  of  ~0.5  nm.  This  GB  solidus  line,  where  presumably  quasi-liquid  IGFs  start  to  form,  is 
marked  as  a  thick  red  dashed  line  in  Fig.  C.  10.  Since  the  densification  represents  an  average 
effect  of  many  GBs  with  different  crystallographic  misorientations  and  inclinations,  a  gradual 
transition  in  GB  diffusivity  was  expected  and  observed.  In  contrast,  the  measured  (effective)  GB 
diffusivities  only  jump  moderately  at  the  bulk  solidus  line  where  the  bulk  liquid  appears;  these 
are  labeled  by  the  blue  stars  in  Fig.  C.  10  (A)  and  (E).  Furthermore,  below  the  bulk  peritectic 
temperature,  there  is  no  discontinuity  in  the  measured  GB  diffusivity  at  the  bulk  solvus  line, 
which  is  labeled  by  the  light  purple  star  in  Fig.  C.  10(B).  This  result  is  fully  expected. 


Bulk  Diffusion 


Figure  C.10  All  trends  in  the  measured  GB  diffusivities  can  be  explained  based  on  the  computed 

GB  diagram  for  Ni-doped  Mo. 

The  measured  GB  diffusivity  increases  with  increasing  temperature  for  pure  Mo,  as  shown  in 
Fig.  C.  10(C).  This  trend  is  well  expected.  Flowever,  the  measured  GB  diffusivity  appears  to 
decrease  with  increasing  temperature  at  a  constant  Ni  fraction  of  0.5  %  (which  lies  in  the  single¬ 
phase  region  for  all  temperatures),  as  shown  in  Fig.  C.  10(D).  While  this  result  is  highly  counter¬ 
intuitive,  it  can  be  fully  explained  from  the  computed  GB  diagram.  This  abnormal  behavior  of 
decreasing  GB  diffusivity  with  increasing  temperature  is  due  to  the  fact  that  Ni  has  a  “retrograde” 
solubility  in  the  Mo  based  BCC  phase.  With  increasing  temperature,  this  fixed  composition  of 
0.5  %  moves  away  to  the  metastable  extension  to  the  solidus  line  (i.e.,  the  black  dotted  line  in  the 
computed  GB  diagram  in  Fig.  C.10).  Consequently,  the  free  energy  penalty  for  forming  an 
undercooled  liquid  increases  with  increasing  temperature  (for  fixed  0.5  %  Ni).  Consistently,  the 
computed  TL  decreases  with  increasing  temperature;  this  is  shown  in  the  inset  of  Fig.  C.  1 0(D).  In 
conclusion,  a  GB  can  “solidify”  with  increasing  temperature!  The  above  comparisons  support  the 
correctness  and  usefulness  of  our  thermodynamic  models  and  the  computed  GB  diagram. 
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(C.5)  Tables  of  Densification  Data 


(A)  Densification  data  of  specimens  that  were  isothermally  sintered  at  1495  °C. 


Linear  shrinkage,  AL/L0 

Time,  mins 

0 

10 

30 

60 

120 

240 

Pure  Mo 

0.0031 

0.0105 

0.0204 

0.0307 

0.0393 

0.0464 

0.25%Ni 

0.0126 

0.0204 

0.0235 

0.0272 

0.0314 

0.0377 

0.50%Ni 

0.0000 

0.0199 

0.0288 

0.0320 

0.0377 

0.0477 

0.75%Ni 

0.0047 

0.0372 

0.0419 

0.0435 

0.0451 

0.0514 

1.00%Ni 

0.0078 

0.0351 

0.0413 

0.0466 

0.0529 

0.0633 

1.50%Ni 

0.0031 

0.0414 

0.0503 

0.0582 

0.0687 

0.0786 

1.75%Ni 

0.0031 

0.0713 

0.0781 

0.0891 

0.0959 

0.1101 

2.00%Ni 

0.0047 

0.0801 

0.0926 

0.1091 

0.1179 

4.00%Ni 

0.0031 

0.1148 

0.1234 

0.1344 

0.1384 

Relative  density,  p 

Pure  Mo 

0.620 

0.631 

0.643 

0.649 

0.660 

0.666 

0.25%Ni 

0.616 

0.629 

0.634 

0.641 

0.649 

0.661 

0.50%Ni 

0.588 

0.621 

0.638 

0.645 

0.656 

0.674 

0.75%Ni 

0.615 

0.673 

0.682 

0.685 

0.688 

0.702 

1.00%Ni 

0.613 

0.661 

0.674 

0.684 

0.696 

0.717 

1.50%Ni 

0.599 

0.667 

0.683 

0.696 

0.716 

0.738 

1.75%Ni 

0.605 

0.727 

0.747 

0.772 

0.789 

0.819 

2.00%Ni 

0.586 

0.727 

0.751 

0.792 

0.816 

4.00%Ni 

0.581 

0.803 

0.839 

0.862 

0.874 

73 


(B)  Densification  data  of  specimens  that  were  isothermally  sintered  at  1450  °C. 


Linear  shrinkage,  AL/L0 

Time,  mins 

0 

10 

30 

60 

120 

240 

Pure  Mo 

0.0031 

0.0110 

0.0142 

0.0183 

0.0236 

0.0283 

0.25%Ni 

0.0094 

0.0162 

0.0189 

0.0236 

0.0273 

0.0330 

0.375%Ni 

0.0031 

0.0162 

0.0194 

0.0236 

0.0288 

0.0367 

0.50%Ni 

0.0094 

0.0215 

0.0241 

0.0283 

0.0325 

0.0383 

0.75%Ni 

0.0047 

0.0288 

0.0325 

0.0425 

0.0461 

0.0571 

1.00%Ni 

0.0047 

0.0356 

0.0388 

0.0477 

0.0519 

0.0718 

1.50%Ni 

0.0031 

0.0540 

0.0660 

0.0776 

0.0833 

0.0894 

1.75%Ni 

0.0031 

0.0692 

0.0854 

0.1006 

0.1127 

0.1242 

2.00%Ni 

0.0031 

0.0655 

0.0833 

0.0980 

0.1111 

0.1237 

4.00%Ni 

0.0047 

0.0949 

0.1127 

0.1263 

0.1357 

0.1426 

Relative  density,  p 

Pure  Mo 

0.616 

0.630 

0.636 

0.643 

0.653 

0.661 

0.25%Ni 

0.613 

0.622 

0.627 

0.634 

0.642 

0.652 

0.375%Ni 

0.602 

0.624 

0.629 

0.636 

0.646 

0.660 

0.50%Ni 

0.625 

0.642 

0.645 

0.652 

0.660 

0.672 

0.75%Ni 

0.609 

0.650 

0.657 

0.673 

0.683 

0.706 

1.00%Ni 

0.623 

0.664 

0.668 

0.687 

0.694 

0.736 

1.50%Ni 

0.601 

0.698 

0.720 

0.746 

0.757 

0.777 

1.75%Ni 

0.607 

0.732 

0.767 

0.803 

0.834 

0.868 

2.00%Ni 

0.605 

0.724 

0.762 

0.798 

0.830 

0.865 

4.00%Ni 

0.606 

0.782 

0.828 

0.866 

0.892 

0.911 
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(C)  Densification  data  of  specimens  that  were  isothermally  sintered  at  1400  °C. 


Linear  shrinkage,  AL/L0 

Time,  mins 

0 

10 

30 

60 

120 

240 

480 

960 

Pure  Mo 

0.0016 

0.0094 

0.0142 

0.0173 

0.0200 

0.0241 

0.0294 

0.0346 

0.25%Ni 

0.0079 

0.0157 

0.0189 

0.0236 

0.0272 

0.0314 

0.0377 

0.0451 

0.375%Ni 

0.0031 

0.0157 

0.0200 

0.0252 

0.0303 

0.0357 

0.0404 

0.0477 

0.50%Ni 

0.0079 

0.0220 

0.0241 

0.0309 

0.0362 

0.0419 

0.0482 

0.0556 

0.75%Ni 

0.0047 

0.0299 

0.0330 

0.0362 

0.0430 

0.0519 

0.0618 

0.0781 

1.00%Ni 

0.0031 

0.0383 

0.0461 

0.0608 

0.0645 

0.0708 

1.50%Ni 

0.0031 

0.0550 

0.0671 

0.0833 

0.0980 

0.1132 

1.75%Ni 

0.0031 

0.0671 

0.0875 

0.1001 

0.1132 

0.1268 

2.00%Ni 

0.0031 

0.0629 

0.0797 

0.0954 

0.1101 

0.1242 

4.00%Ni 

0.0031 

0.0928 

0.1174 

0.1316 

0.1405 

0.1520 

Relative  density,  p 

Pure  Mo 

0.618 

0.632 

0.641 

0.647 

0.651 

0.659 

0.668 

0.679 

0.25%Ni 

0.619 

0.630 

0.635 

0.644 

0.650 

0.657 

0.670 

0.684 

0.375%Ni 

0.607 

0.628 

0.637 

0.646 

0.656 

0.665 

0.675 

0.689 

0.50%Ni 

0.624 

0.645 

0.648 

0.659 

0.669 

0.680 

0.693 

0.709 

0.75%Ni 

0.623 

0.656 

0.661 

0.667 

0.680 

0.699 

0.720 

0.754 

1.00%Ni 

0.621 

0.671 

0.682 

0.712 

0.721 

0.733 

1.50%Ni 

0.602 

0.699 

0.723 

0.759 

0.794 

0.834 

1.75%Ni 

0.601 

0.723 

0.765 

0.799 

0.830 

0.869 

2.00%Ni 

0.606 

0.718 

0.754 

0.790 

0.827 

0.866 

4.00%Ni 

0.577 

0.743 

0.804 

0.841 

0.869 

0.904 
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(D)  Densification  data  of  specimens  that  were  isothermally  sintered  at  1350  °C. 


Linear  shrinkage,  AL/L0 

Time,  mins 

0 

20 

60 

120 

240 

480 

940 

Pure  Mo 

0.0047 

0.0083 

0.0132 

0.0168 

0.0220 

0.0247 

0.0288 

0.25%Ni 

0.0110 

0.0162 

0.0199 

0.0225 

0.0278 

0.0335 

0.0383 

0.375%Ni 

0.0031 

0.0157 

0.0199 

0.0241 

0.0320 

0.0404 

0.0451 

0.50%Ni 

0.0079 

0.0220 

0.0278 

0.0341 

0.0425 

0.0503 

0.0540 

0.75%Ni 

0.0063 

0.0262 

0.0419 

0.0618 

0.0734 

0.0802 

1.00%Ni 

0.0031 

0.0294 

0.0440 

0.0681 

0.0839 

0.0959 

1.25%Ni 

0.0063 

0.0288 

0.0451 

0.0723 

0.0922 

0.1085 

2.00%Ni 

0.0031 

0.0330 

0.0482 

0.0749 

0.0970 

0.1116 

4.00%Ni 

0.0031 

0.0346 

0.0561 

0.0828 

0.1053 

0.1216 

Relative  density,  p 

Pure  Mo 

0.611 

0.617 

0.625 

0.632 

0.641 

0.647 

0.655 

0.25%Ni 

0.627 

0.633 

0.640 

0.645 

0.655 

0.666 

0.675 

0.375%Ni 
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0.633 

0.641 

0.649 

0.664 

0.680 

0.690 

0.50%Ni 

0.631 

0.650 

0.658 

0.671 

0.687 

0.702 

0.710 

0.75%Ni 

0.612 

0.646 

0.677 

0.717 

0.743 

0.758 

1.00%Ni 

0.603 

0.650 

0.678 

0.729 

0.766 

0.792 

1.25%Ni 

0.612 

0.651 

0.683 

0.739 

0.787 

0.826 

2.00%Ni 

0.612 

0.666 

0.694 

0.753 

0.804 

0.842 

4.00%Ni 

0.603 

0.661 

0.702 

0.761 

0.815 

0.854 
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(E)  Densification  data  of  specimens  that  were  isothermally  sintered  at  1300  °C. 


Linear  shrinkage,  AL/L0 
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0.0215 
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0.0922 

1.75%Ni 
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0.643 

0.650 

0.660 

0.669 
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0.631 
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0.643 

0.652 

0.664 

0.674 

0.50%Ni 

0.619 

0.640 
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0.673 

0.696 

0.708 

0.717 

0.75%Ni 

0.606 

0.637 

0.668 

0.699 

0.743 

0.782 

1.75%Ni 

0.601 

0.642 

0.672 

0.700 

0.749 

0.793 

4.00%Ni 

0.584 

0.620 
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(F)  Densification  data  of  specimens  that  were  isothermally  sintered  at  1220  °C. 


Linear  shrinkage,  AL/L0 
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0.0110 
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0.0346 
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0.0354 
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0.0634 

0.0818 

2.00%Ni 
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0.0236 
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0.0542 
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0.739 

0.783 

1.00%Ni 

0.598 
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0.649 

0.661 

0.704 

0.743 

2.00%Ni 
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0.625 

0.656 

0.684 

0.720 
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(G)  Densification  data  of  specimens  that  were  isothermally  sintered  at  1160  °C. 


Linear  shrinkage,  AL/L0 
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(H)  Densification  data  of  specimens  that  were  isothermally  sintered  at  1 100°C. 
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